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Abstract 

We explore various aspects of the motion of heavy quarks in strongly-coupled 
gauge theories, employing the AdS/CFT correspondence. Building on earlier 
work by Mikhailov, we study the dispersion relation and energy loss of an 
accelerating finite-mass quark in N = 4 super- Yang-Mills, both in vacuum and 
in the presence of a thermal plasma. In the former case, we notice that the 
application of an external force modifies the dispersion relation. In the latter 
case, we find in particular that when a static heavy quark is accelerated by an 
external force, its rate of energy loss is initially insensitive to the plasma, and 
there is a delay before this rate approaches the value derived previously from 
the analysis of stationary or late-time configurations. 

Following up on work by Herzog et al, we also consider the evolution of a 
quark and antiquark as they separate from one another after formation, learn- 
ing how the AdS/CFT setup distinguishes between the singlet and adjoint 
configurations, and locating the transition to the stage where the deceleration 
of each particle is properly accounted for by a constant friction coefficient. Ad- 
ditionally, we examine the way in which the energy of a quark-antiquark pair 
moving jointly through the plasma scales with the quark mass. We find that 
the velocity-dependence of the screening length is drastically modified in the 
ultra-relativistic region, and is comparable with that of the transition distance 
mentioned above. 
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1 Introduction and Summary 
1.1 Brief overview of earlier work 

In the last couple of years, an intense research effort has been directed toward the 
use of the AdS/CFT correspondence [HEIE] to study parton energy loss in strongly- 
coupled thermal plasmas. This endeavor is motivated mostly by the quest to under- 
stand the quark-gluon plasma (QGP) [I] produced at RHIC [5] and LHC [6], and 
was stimulated by the pioneering works [3 El US HQ], which were in turn encouraged 
by the success of the earlier viscosity cakulation7[IH [12] The phenomenological 
literature on energy loss is enormous; for reviews, see, e.g., [T4l [T5] . 

The drag force experienced by a heavy quark traversing an M = 4 super- Yang- 
Mills (SYM) plasma was determined in E], via its dual description as a string 
moving on an AdS-Schwarzschild background. The closely related diffusion coefficient 

1 Recently it has been argued [13] that, at least at weak coupling, there in fact exists a direct 
link between viscosity and jet quenching. 
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was obtained independently in [9]. These seminal papers have been generalized and 
elaborated on in a vast number of posterior contributions rT6 |[T7l[T5|[T^l[2m[2Tl[2^| [23]. 
including in particular comparisons with the corresponding weakly-coupled results 
[23], as well as extensive analyses of the energy- momentum tensor which paint a 
detailed and beautiful picture of the directionality of energy flow away from the 
moving quark [25] . 

A second line of development originated from the work [TU], whose authors pro- 
posed a recipe for the jet-quenching parameter q used in some phenomenological 
models of energy loss [H], and employed it in the context of M = 4 SYM. Their 
calculation has been extended in a number of directions in [261 El 1211, but also 
questioned in various ways in [28l [181 [T9l [291 [201 [30] . 

A third line was initiated in [321 [28], which studied mesons moving through an 
M = 4 SYM plasma and obtained the corresponding quark-antiquark potential and 
screening length, using the dual portrayal in terms of a string that moves on an AdS- 
Schwarzschild background and has both of its endpoints on the boundary. As was 
emphasized in [32] , the outcome would be expected to have implications for the issue 
of quarkonium suppression in the QGP. Related results were obtained independently 
in [21], which determined the spectrum of spinning mesons in the confining chiral 
gauge theory dual to the Sakai-Sugimoto model [33]. Various interesting extensions 
and refinements of these calculations have been reported in [331 EH [361 EH EH [211 EH 
[221 SO]. 

A notable feature is that, in contrast with the quark probes considered in [3, [HI 
El [TO] , and the gluon probes studied in [HI [23] , mesons do not feel a drag force as 
they move through the plasma [2H EH [2S] • This is because they are color- neutral, 
and therefore incapable of setting up the long-range gluonic field profiles that could 
transport energy away from themo Indeed, in [T8l [22] it has been shown that the 
other obvious color-neutral probe, the baryon, likewise experiences no drag. 

Besides the transport and jet quenching properties of the plasma, there have also 
been interesting studies of photoemissivity [33] and deep inelastic scattering [441 140]. 
Some of the topics we have briefly enumerated here have been reviewed in more depth 
in [33 [36]. 

1.2 Motivation and main results 

Naturally, the initial papers [Zl [HI El Ell E21 [28] carried out their calculations under a 
number of simplifying assumptions. First and foremost among these is of course the 
use of an M = 4 SYM plasma as a toy model for the real-world QGP. A number of 
extensions to other theories that are in certain ways closer to QCD have already been 
cited above, and a few others can be found in [37] . The issue of how best to compare 
the M = 4 SYM and QCD parameters was discussed in [31] . Also important is the 
restriction to an infinite static plasma. Efforts to examine the case where the plasma 

2 The color field profile set up by a meson has been explicitly determined at zero temperature in 

[H1H2]. 
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is expanding and/or has a finite spatial extent have been made in [4*9] . 

The plasma produced at RHIC or LHC has in addition a finite temporal extent, so 
another issue that could be significant is the limitation of the energy loss computations 
[3 El El [10] to the stationary or late-time regime. The actual drag coefficient might 
be expected to differ from the value obtained in these works in a situation where the 
quark moving through the plasma is accelerating, or in the initial period following its 
production within the thermal medium. This point has been emphasized from the 
phenomenological perspective (in the context of collisional energy loss) in [5U], where 
it was argued that, after a quark-antiquark pair is produced, there exists a significant 
delay before the rate of energy dissipation coincides with the result relevant for the 
stationary case. Later works [51] have called into question the actual duration, but 
not the existence, of this retardation effect. 

The estimates in \5U\ 1ST] are based on perturbative calculations, so it is interesting 
to inquire into this effect in the strongly-coupled systems available to us through the 
AdS / CFT correspondence. It was this question that got us started on the investiga- 
tion that led to the present work, which over time has expanded toward various other 
fronts. For simplicity, we have carried out our calculations in the context of an M = 4 
plasma, but we expect most of our qualitative conclusions to apply more generally. 

Now, of course, the restriction in [8] to the stationary and asymptotic cases 
was not made gratuitously, but was necessary in order to gain analytic control on the 
problem of energy loss. Away from these regimes, the evolution of the quark in the 
thermal plasma, or equivalently, of the string on the AdS-Schwarzschild geometry, 
is complicated, and one must resort to a numerical analysis (aspects of which were 
explored already in [7]). As we will examine closely in due course, the interpretation 
of the outcome of such an analysis is encumbered by several difficulties, chief among 
which is our ignorance of the thermal dispersion relation for the quark. 

In Section [2] we will gain some perspective on these issues by turning off the tem- 
perature T, and analyzing first the evolution of an isolated quark in vacuum, which 
is of interest in its own right, and easier to interpret because in this context the form 
of the dispersion relation is fixed by Lorentz invariance. We begin in Section 12.11 by 
reviewing a remarkable paper by Mikhailov [52], who constructed an analytic embed- 
ding for the string dual to an infinitely-massive quark in M = 4 SYM that follows 
an arbitrary timelike trajectory, and extracted from it a rate of energy loss, Eq. 0, 
which turned out to agree with the Lienard formula from classical electrodynamics! 
In the derivation of his result, Mihailov disregarded a total derivative, which we show 
in ([TO]) to give precisely the expected dispersion relation for the quark. From the re- 
sults of [52], then, we are able to draw two important lessons: first, that, at any given 
time, the energy of the string includes not only the portion intrinsic to the quark, but 
also the part that has been radiated by the quark throughout its previous history; 
second, that the quark's dispersion relation arises from a total derivative that ends 
up being evaluated at the string endpoint. 

It is natural to wonder how Mikhailov's results generalize to the case where the 
quark has a finite mass, which requires the introduction of probe D7-branes on which 
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the string can end [53] . We examine this issue (still at T = 0) in Section 12.21 Inter- 
estingly, we discover that the resulting quark dispersion relation, Eq. (1T91 , as well as 
the rate of energy loss, Eq. (ITS]) , depend on the external force F exerted on the quark, 
or equivalently, on the string embedding parameter X', which in the gauge theory 
controls the shape of the 'gluon cloud' surrounding the quark. The dependence is 
such that the energy and momentum of the quark reduce to the familiar expressions 
when F — > 0, and on the other hand diverge as the force approaches its critical value 
(i.e., the value beyond which F would be strong enough to nucleate quark- ant iquark 
pairs out of the vacuum). We finalize our zero-temperature analysis in Section [2~3| 
noting that the process of energy loss is accompanied by the formation of an event 
horizon (and a stationary limit curve) on the string worldsheet, as depicted in Fig. [TJ 
Armed with the intuition afforded to us by Mikhailov's construction, we proceed 
in Section [3] to the finite-temperature case. In this setting, we are of course limited 
by the fact that the general solution to the string equation of motion is not known 
analytically. Nonetheless, we argue that Mikhailov's method should admit a T > 
generalization, and are able to show in Section 13.11 that this is indeed true for the 
only thermal solution that is thus far available in closed form, which corresponds to 
the stationary quark configuration studied in [7J, [S] . The resulting dispersion relation, 
seen in the second and third line of Eq. (151]) ; contains a novel feature that we argue to 
apply for all finite-temperature configurations, including the quark at rest, Eq. (138]) : 
it receives a contribution not only from the string endpoint located on the D7-branes, 
which is directly dual to the quark, but also from the endpoint located at the black 
hole horizon, which, as we explain, encodes the initial conditions for the joint quark 
+ plasma system. 

In Section 13.21 we consider the more general case where the quark accelerates 
within the plasma. After reviewing the work done in this context by the authors of 
[TJ, we pick up precisely where they left off, integrating the string equation of motion 
numerically for an initially static quark that is accelerated by an external force over a 
finite period of time and is thereafter released. As shown in Figs. I2~f51ffl we find that 
under such conditions, and for values of the mass in the neighborhood of the charm 
quark, there exists a period after release where the quark dissipates energy at a rate 
that is substantially smaller than the stationary/asymptotic result (1571) obtained in 
[7JE]. Additionally, as seen in Fig. El the rate of energy loss in the initial stage where 
the quark is pushed by the external force can be almost completely accounted for by 
the generalized Lienard formula (1 18]) . describing radiation in vacuum. In (152]) and 
the paragraphs immediately following, we discuss the form of the thermal dispersion 
relation for the quark subject to the above initial conditions, incorporating the in- 
dependence expected from our zero-temperature results, and explaining in detail the 
relation to the relativistic expression (1151) proposed in \7\. In Section [3731 we then carry 
over the discussion of the worldsheet black hole from Section [2731 to the T > context, 
obtaining a time-dependent analog of the Schwarzschild black hole encountered for 
the stationary case in [191 120], as schematized in Fig. [7J 

Having explored in some depth the effect of the acceleration on the rate of energy 
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loss for an isolated quark, in Section[4]we introduce the other element of realism whose 
importance was highlighted by the phenomenological studies [SU] EEL an d examine a 
quark that is produced together with its corresponding antiquark at some finite time 
within the thermal medium. The gravity description of this system involves a string 
with both of its endpoints on the D7-branes, at initially coincident positions. As we 
review in Section I4.1[ an initial exploration of this type of configuration was carried 
out in [7|. In Section H~2l we show how the one-parameter family of initial conditions 
considered in that work can be generalized to describe situations with vastly different 
patterns of excitations for the initial gluonic fields, and more importantly, to the case 
where the newly- formed quark and antiquark transform in the adjoint instead of the 
singlet representation of the color gauge group. Curiously, we find that it is only in 
the adjoint case that the initial quark velocity v is freely adjustable. For the singlet 
case, as in [7] one must necessarily have v = v m , where v m is the velocity (I58|) that 
appeared in different manifestations in the works [3~T] 1321 128] [T9l [20] . and on the string 
theory side of the duality corresponds to the proper velocity of light at the location 
of the string endpoint [35]. This identification strongly suggests that v m should be a 
limiting velocity, as is implicit in at least [22 EH ED] , and is discussed more explicitly 
in Section 14.31 (as well as in the very recent works [3"5l [30] , which appeared while this 
paper was in preparation). 

In Section H~4l we study the transition of the q-q pair to the asymptotic regime de- 
scribed by a constant friction coefficient, which as shown in Fig. [10] applies uniformly 
to all different types of initial conditions. For singlet configurations, the initial evo- 
lution of the quark is of course drastically affected by the presence of the antiquark. 
In Fig. [TT] we determine the velocity dependence of the transition distance beyond 
which the quark is effectively in the asymptotic regime. For adjoint configurations, 
on the other hand, the interactions between the quark and antiquark are suppressed 
at large N, and so the two members of the pair evolve independently from the start. 
In this case we find that the transition distance is essentially zero. 

A natural question is whether the transition to the asymptotic regime occurs 
right after the quark and antiquark are screened from one other by the plasma. This 
requires a determination of the corresponding quark-antiquark potential, and more 
specifically the screening length, a problem that we turn to in Section [5j After 
reviewing and comparing in Section 15.11 the results obtained for infinitely-massive 
quarks in [32] EE], we generalize to the case of finite mass, first at zero temperature 
in Section 15.21 and then at finite temperature in Section 15.31 The resulting potentials 
are shown in Figs. [T3l [T4l As expressed in f[84l . ([88]) . (1941 . they are found to be linear 
instead of divergent when the quark and antiquark approach one another, signaling 
the fact that the color sources in this case are no longer pointlike. 

In the finite temperature case, the potential implies the velocity-dependence of 
the screening length L max presented in Figs. [T5l[T6l As seen there, compared to the 
infinitely- massive case examined in [32], [28], there is a drastic modification of the 
behavior at high velocities, which are now bounded by v m instead of 1. The in- 
dependence near this limit can be determined analytically and takes the form ([97]) . 
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instead of the formula (180|) obtained in [32J. Over the whole range < v < v m , and for 
masses similar to the charm quark, the behavior can be relatively well approximated 
by dHS]), which is the obvious generalization of the fit (ITS]) proposed in [28] ■ 

We end the paper by comparing in Section 15.41 the screening length against the 
transition distance determined in 14.41 For singlet configurations, we find that the mag- 
nitude and velocity-dependence of these two separations are comparable, as shown 
in Fig. [TTJ Notice that this is in spite of the fact that the two relevant string con- 
figurations are rather different. A similar statement can be made for adjoint config- 
urations, where both the screening and transition lengths are essentially zero. We 
conclude then that, for q-q pairs created within the plasma, the transition to the 
asymptotic, constant-drag-coefficient regime takes place immediately after the quark 
and antiquark lose contact with one another. That is to say, there is no intermediate 
stage where the quark and antiquark decelerate independently from one another at a 
rate that differs substantially from the asymptotic result of [8J . 

Our analysis in Section 13.21 demonstrates that the gluonic fields around a quark 
can be disturbed by the application of an external force, to the point of producing a 
significant modification of the energy dissipation rate in the period immediately fol- 
lowing release. By analogy with the results we obtained for singlet quark- antiquark 
configurations in Section 14. 4[ we interpret this to mean that the quark has to escape 
far enough from the disturbed region in order to be screened from its effect. In any 
event, given that the actual experimental situation resembles the setup of Section [4] 
much more closely than that of Section 13.21 the main overall lesson for QGP phe- 
nomenology would appear to be that, beyond an initial period which is controled by 
the screening length, and where the evolution can be modeled relatively well as taking 
place in vacuum (and, as such, would be present also in proton-proton collisions), the 
stationary/asymptotic rate of energy dissipation determined in [7], [8J, [9] gives a good 
approximation to the actual time-dependent dynamics. 

It is worth emphasizing that the inferences made in this paper regarding disper- 
sion relations and energy loss rates are based on the natural split achieved in [52] of 
the total energy of the string. The latter is conserved on the fixed AdS-Schwarzschild 
background, but would of course decrease steadily if we take into account the grav- 
itational (and dilatonic, etc.) radiation given off by the string in the course of its 
evolution. Through the GKPW recipe for correlation functions [2], it is this radiation 
(or, more precisely, the full metric perturbation produced by the string), evaluated at 
the AdS boundary, that determines the expectation value of the gauge theory energy- 
momentum tensor. As has been meticulously studied in [25] , this tensor contains not 
just the gross information about the total energy loss rate, but even the fine details 
about the directionality of the flow and the relative weight of the various dissipation 
channels. It would therefore be very interesting to examine more closely this link 
between dissipated energy as encoded on the string and on the gravitational field 
ultimately generated by it. 
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2 Single Quark Evolution: Zero Temperature 



To analyze the motion of heavy quarks in a strongly-coupled M = 4 SU (N c ) SYM 
plasma with coupling gyu and temperature T, one must follow the evolution of open 
strings that end on a stack of Nf D7-branes [53] living on the (AdS-Schwarzschild)s x 
S 5 geometry 

ds 2 = G mn dx m dx n = (-hdt 2 + dx 2 + + R 2 dn 5 , (1) 

z l \ h J 

z A i? 4 1 

h = I - — , -j±= gy M N c = A , z h = — , 

Z h l s 711 

where l s denotes the string length. In our work the string will be taken to lie at a 
fixed position on the S 5 (consistent with the corresponding equations of motion), so 
the angular components of the metric will not play any role. The D7-branes cover the 
four gauge theory directions t, x, and extend along the radial AdS direction up from 
the boundary at z = to a position where they 'end' (meaning that the S 3 C S 5 that 
they are wrapped on shrinks down to zero size), whose location z = z m is related to 
the quark mass in a way that we will specify below. 

From the gauge theory perspective, the introduction of the D7-branes in the back- 
ground (0Q) is equivalent to the addition of Nf hypermultiplets in the fundamental 
representation of the SU(N C ) gauge group, breaking the supersymmetry down to 
M = 2. These are the degrees of freedom that we refer to as 'quarks,' even though 
they include both spin 1/2 and spin fields. For Nf <C N c , the backreaction of the 
D7-branes on the geometry can be sensibly neglected; in the field theory this corre- 
sponds to working in a 'quenched' approximation which disregards quark loops (as 
well as the positive beta function they would generate). 

The string dynamics follows as usual from the Nambu-Goto action 

Sng = ~7^7 / d2(X V~ det 9ab = 7j^7 / dV£ NG , (2) 

where g a b = d a X m dbX n G mn (X) (a,b = 0,1) denotes the induced metric on the 
worldsheet. In the static gauge a = t, a 1 = z, and for motion and deformation of 
the string purely along direction x = x 1 , the non-zero canonical momentum densities 
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11° = dC N G/d(d a X^) are given by 



hx' 2 + 1 



K = . , (3) 



z 2 y 1 + /zA 


X 2 






z 2 hsJl + hX' 2 - 








z 2 y / l + /iX' 2 - 


) 

X 2 

h 


/iX' 





z\ll + hX' 2 -^ 



where of course X = d t X(t, z), X' = d z X(t, z). Notice that, due to our normalization 
of £ng> the 11° must be multiplied by R 2 /2ira' = \f~\/2n to obtain the physical energy 
and momentum densities. 

In the present section we will restrict attention to the case of vanishing temper- 
ature (zh — > oo), in which case we are left in flJJ) with a pure AdS geometry, and 
the D7-brane parameter z m is inversely proportional to the Lagrangian mass of the 
quark, 

z m = -z • 4 

Zirm 

A quark that accelerates in vacuum would be expected to emit chromoelectro- 
magnetic radiation. This problem has been examined from the classical perspective 
in [51], and quantum-mechanically at weak coupling in, e.g., [55]. First steps towards 
a strong-coupling analysis by means of the AdS/CFT correspondence were taken in 
[JT] , which employed tools developed in [56] to study the dilatonic waves given off by 
small fluctuations on a radial string in AdSs, and infer from them the profile of the 
gluonic field (TrF 2 (x)} in the presence of a quark undergoing small oscillations. The 
results of [JT] painted an interesting picture of the propagation of nonlinear waves in 
M = 4 SYM, but did not allow a definite identification of waves with the l/\x\ falloff 
associated with radiation. Very recently, this falloff has been successfully detected in 
the same setup as [JT] through a calculation of the energy-momentum tensor (T^ u ) 
|57j . which appeared while this paper was in preparation. 



2.1 Infinite mass 

The first definite characterization of the radiation rate off an accelerating quark was 
found much later than [JT], and by a completely different route, in an important 
paper by Mikhailov [52] . Remarkably, this author was able to solve the full nonlinear 
equation of motion for a string on AdSs, for an arbitrary timelike trajectory of the 
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string endpoint dual to a heavy quark! In terms of the coordinates used in ([Q) (where 
for now h = 1), his solution is 



X»(t,z) = z^^ + x»(t) , (5) 

with [j, = 0,1, 2, 3, and x M (r) the worldline of the string endpoint at the AdS boundary— 
or, equivalently, the worldline of the dual, infinitely massive, quark — parametrized 
by the proper time r defined through r) fll/ x ll x 1/ = — 1, where x M = dx^/dr. Equation 
([5]) displays the string worldsheet as a ruled surface in AdSs, spanned by the straight 
lines at constant r. 

Combining ([T]) and §5§, the induced metric on the worldsheet is found to be 

d2 e>2 

9rr = —{z x - 1), g zz = 0, g ZT = -, 

implying in particular that the constant-r lines are null, a fact that plays an important 
role in Mikhailov's construction. In the solution ([5]), the behavior at time t = X°(t, z) 
of the string segment located at radial position z is completely determined by the 
behavior of the string endpoint at a retarded time t ret (t, z) obtained by projecting back 
toward the boundary along the null line at fixed r. From the /x = component of 
parametrizing the quark worldline by x°(t) instead of r, and using dr = \/l — v 2 dx°, 
where v = dx/dx°, this amounts to 

* = z n 1 ^9 + trct ' ( 6 ) 
VI — v 

where the endpoint velocity v is meant to be evaluated at t vet . In these same terms, 
the spatial components of ((SJ) can be formulated as 

-> v 

X(t, z) = z + x(t ret ) = {t- t iet )v + x(t iet ) . (7) 

y/l — V z 

Working in the static gauge a = t, a 1 = z, the total energy of a string that 
extends all the way down to the boundary — i.e., with z m = 0, corresponding to an 
infinitely massive quark — follows from the Nambu-Goto action as 



E(t) 



Vx f°° dz 1 + (Hr) 



271 Jo z2 li _ (dxY + (ox) 2 _ (dx) 2 (dx) 2 + (ox . dx\ 2 



at J ^ yaz j yat j \dz j ^ yat dz J 



Using ([6]) and ([7]), Mikhailov was able to reexpress this energy (via a change of 
integration variable z —>■ t TCt ) as a local functional of the quark trajectory, 

\A a 2 - \v x a] 2 



E(t) = — / dt Iet 1 J + E q (v(t)) , (9) 
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where of course a = dv/dx°. The second term in the above equation arises from a 
total derivative that was not explicitly written down by Mikhailov, but can easily be 
worked out to be 

E 9 W = %-(-==-) =7m f (10) 

which gives the expected Lorentz-invariant dispersion relation for the quark. The en- 
ergy split achieved in therefore admits a clear and pleasant physical interpretation: 
E q is the intrinsic energy of the quark at time t, and the integral over t ret encodes the 
accumulated energy lost by the quark over all times prior to t. No less remarkable is 
the fact that the rate of energy loss for the quark in this strongly-coupled non-Abelian 
theory is found to be in precise agreement with the standard Lienard formula from 
classical electrodynamics! 

For the momentum of the string, Mikhailov analogously found 

V~\ /** ,. a 2 — \v x a} 2 



where the second term can be worked out to be 

_ _ VA / v 1 



jmv , (12) 



and encodes the momentum intrinsic to the quark. The last two equations also follow 
from and ( flOl) through Lorentz invariance. 



2.2 Finite mass 

It is interesting to consider how these results are modified in the case z m > 0, where 
the mass m of the quark given by (j4j) is large but not infinite. As we have reviewed 
above, in Mikhailov's original solution (jSJ) the evolution of the string at any radial 
position z follows from knowledge of the trajectory of the endpoint at the AdS bound- 
ary. For a finite-mass quark, we ought to impose boundary conditions on the string 
not at z = but at z = z m : given the worldline x(t) of the quark, we must require 
that the string worldsheet satisfy X(t, z m ) = x(t). Moreover, we need only determine 
the behavior of the string in the region z > z m . 

The required physical solution can of course be viewed as merely the z > z m 
portion of one particular instance of the general solution found by Mikhailov. Our 
task is therefore to reexpress (|SD in terms of the data x(t) at the new boundary 
z = z m , instead of the (now merely auxiliary) data at the AdS boundary z = 0, 
which we will henceforth distinguish with a tilde: x(t). For simplicity, we will carry 
out this translation explicitly only in a setup where the quark moves purely along 
direction x = x 1 , which is all that we will need for our analysis in subsequent sections. 
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It follows from and (J7|) that, at any given point (t, z) on the string worldsheet, 



, dz , 

at = , + dt mt 



vT 



r,2 



Ct^L = - -|- Utrpf 



1_ £2)3/2 ' 



(13) 



+ 7^ TTTT^T +^ 



[ V /T^2 (1-^)3/2 



where u, a denote the velocity and acceleration at the point (t = t rct ,z = 0) on the 
AdS boundary obtained by projecting back from (t, z) along a null trajectory. From 
(113"]) we can deduce that 

fdX_\ dz + v(l-v 2 )^ 2 

V dt ) z ~ vaz + (1 - C 2 ) 3 /2 • ( > 

Evaluated at the new boundary z = z m , this formula relates the velocity v = dx/dt = 
dtX(t, z m ) of the actual string endpoint — i.e., the velocity of the finite-mass quark — 
to the velocity v and acceleration a of the 'auxiliary endpoint' at z = 0. Equation 
( TT4l) implies that the quark acceleration a = d 2 x/dt 2 = d 2 X(t, z m ) depends not only 
on v and a, but also on the second time derivative of v. Because of this, it is not 
possible to solve for v and a, the quantities that appear directly in Mikhailov's energy 
formula (Q, in terms of v and a, the data that we would naively expect to suffice to 
characterize the rate of energy loss of the heavy quark. 

On the other hand, from and (J7D we can infer as well that 



fdx\ VT^Pdz 

V dz J t vdz + (1 - v 2 f' 2 ' 1 ' 

so it is certainly possible to solve for v, a in terms of v and d z X(t, z m ). Notice that 
when z m —* 0, we automatically have d z X(t, z m ) — > 0, which explains why this 
parameter was not needed in the description of the infinitely-massive quark. For 
z m > 0, the value of d z X(t, z m ) encodes how much the string tip tilts away from the 
vertical. But what does this mean in gauge-theoretic language? Through a calculation 
of (TrF 2 ) or (T^) in parallel with that of [56l HH [25], the slant of the string will 
have an impact on the shape of the gluonic field profile in the immediate vicinity of 
the heavy quark: for d z X(t,z m ) ^ 0, this profile is not spherically symmetric. We 
find it physically more transparent to express the rate of energy loss in terms not of 
d z X(t, z m ) but of the string momentum density Tl z x given by which by use of (|T4j) 
and (|T5|) can be rewritten as 

K = 7(^)37? • < 16 > 

When evaluated at z = z m , this controls the external force F = (v^A/27r)n^(t, z m ) 
acting on the string endpoint, or equivalently, on the quark. 
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Inverting (fl4|) and (|T6|) . we find 



1 - v n ' 

'l-uWl 



(17) 



r i-zivm 



where we have abbreviated II = H. z x . Using this and ([TBI in (Q, we finally conclude 
that the total energy of the string at time t is given by 



dtzLW 



vUzl 



ztll 2 



As before, the first term represents the accumulated energy lost by the quark at all 
times prior to t: it is the generalization to the m < oo case of the Lienard formula 
([9]) deduced by Mikhailov. The second term again denotes a surface term and gives 
the modified dispersion relation for the finite-mass quark, 



E g (v,F) 



1 - z 2 m vU 



v/A 



27T \ Z y/{\-tf){l-ztfP) 



2vrm 2 - y/XvF 
sjAn 2 m A - a"P 



7m 



(19) 



Starting instead from Mikhailov's formula ffTTl) for the momentum, we find 



dtziU 2 



v — Hz, 



ziU 2 



where 



Pq(v,F) 



v - z 2 Jl 



27T \ Z y/(l-V*)(l-ztfP) 



2nm 2 v - V\F 



7m 



(20) 



(21) 



Notice that dE q /dp q = {2nm 2 v — y/\F) / \2-Km 2 — y/XvF), which for m < 00 and 
F 7^ differs from the result expected for a pointlike quark, dE/dp = v. This reflects 
the fact that the fundamental source dual to a string that terminates at z m > is 
indeed not pointlike. According to the standard UV/IR connection [58], it has a linear 
size of order z m , and it is only because of this extended nature that, as we saw above, 
to characterize its state one needs to specify not only the velocity v but also the 
parameter F (or d z X(t, z m )) that encodes its shape. The crucial point here is that 
the source in question should not be thought of as a bare quark, but as a 'dressed' 
or 'constituent' quark, surrounded by a gluonic cloud with thickness z m [591 [60]. We 
will see more evidence of this in Section 

Another salient feature of the energy and momentum of the quark given by ex- 
pressions f[T9|) and (]2ip is the fact that they both diverge as the value of the external 
force approaches 

2ixm 2 

-fcrit = 7=^ • (22) 

V A 
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The reason for this is easy to understand on the string theory side. To exert a force 
F on the string endpoint, within the D7-branes we must turn on an electric field that 
has strength F Q i = F at z = z m . Working in static gauge, the Born-Infeld Lagrangian 
on the D7-branes is then 



2 



y/-det{g ab + 27ia'F ab ) oc y/-G a G xx - (27ra'F 01 ) 2 = R 2 J ^ - (^ F oi) 



which is real at z = z m only as long as the electric field is below the value F^ = 
\/X/27rz^. Through (gj, this is seen to coincide with the value of the critical force 
The physical origin of this bound is the fact that, for F Q i > F^, the creation of 
open strings is energetically favored, and so the system is unstable. According to (|T9l) 
and (12TI) . then, the energy and momentum of the constituent quark diverge precisely 
at the point where the external force becomes capable of nucleating quark- ant iquark 
pairs. 

2.3 Late-time behavior and worldsheet black hole 

It is interesting to consider the evolution of the string in the case where the quark is 
accelerated by an external force F(t) over some period of time and is then set free 
at a time t re i e ase- If we put F — (II = 0) in (1171) . then a = (the acceleration of 
the auxiliary z = endpoint of the string vanishes), which in turn implies through 
(I14p and ()15p that X — v and X' = at all points on the same null (constant t ret ) 
line, independently of the value of z. In particular, the lower (z = z m ) endpoint of 
the string travels at speed v = X(t, z m ) = v, which means that that, for t > t re i ease , 
the quark moves at constant velocity v, as one would expect given the fact that it is 
in vacuum. 

As time progresses, the null line that departs from the point (t re \ ease , z m ) , which 
according to ([6]) is given by 



^release (^) VI V (t ^release) Zm j (^^) 

reaches further away from the AdS boundary, so there is an increasing portion of the 
string (z m < z < Release (t)) that is completely vertical and moves with the same final 
velocity v. As t — > oo, this vertical segment describes a quark that moves at constant 
speed and is surrounded by a stationary gluonic field profile, related to that of a static 
quark through a Lorentz transformation. 

At any given time t > t rc \ casc , then, the energy and momentum previously radiated 
by the quark are stored as excess energy and momentum on the portion of the string 
above z re i ease (t). This suggests the existence of a geometric region on the worldsheet 
that 'absorbs' the surplus E and p. 

To define this region more precisely, we should note first that, for the spacetime 
metric (OQ) (with h = 1), null curves on the worldsheet obey 

= (24) 
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with the upper or lower sign for the upward- or downward-pointing half of the future 
light cone, respectively. Using (JHJ), (IToT) and (fTTj) . this translates into 




V(i-^ 2 )(i- 



(25) 



and 



(i - n) ^(1-^(1-^) [i - (z* + zl)zlv?\ 



(26) 



zz m U(v - zlU) + (1 - vziTl)^/T^4JP 



where, as before, v and II refer to the velocity and external force at z m , evaluated at 
the retarded time corresponding to the given (t, z). As mentioned earlier, Mikhailov's 
constant-t rct lines (E])-©, and fl23l) in particular, are null, and indeed they can be 
easily seen to satisfy fl25l) . 

If we follow the point z re \ ease (t) as t — > oo (reaching the spacetime horizon for the 
Poincare patch, z — > oo), and then project back along the downward pointing light 
half-cone z^J^t), we delineate a region [z > z^J^t)) on the string worldsheet from 
which, by construction, no signal can escape to the asymptotic region corresponding 
to the final vertical and stationary segment of the string. In other words, the curve 
z rm\\(t) so obtained, which we will henceforth denote by zbh(^), is the event horizon 
of a worldsheet black hole. 

A worldsheet black hole figured prominently in the energy loss analysis of [201 IT9] , 
concerning a quark in a thermal plasma. Those works considered the steady-state 
configuration where the quark moves at a constant velocity v ; as a consequence, the 
black hole they encountered was static, with an event horizon located at the fixed 
position zbh(^) = Zh(l — v 2 ) 1 ^ 4 , precisely the radius that played a crucial role in the 
drag force calculation of [3, IB]- 

In our non-stationary system, on the other hand, the worldsheet black hole is 
dynamical. Let us focus for concreteness on the case where the quark is static up 
to a time t = t gra b, and is then accelerated until t = t re i ease . The location Zbh(^) 
of the event horizon will begin descending from z — > oo even before t = t gra b (in 
anticipation of the disturbance produced by the acceleration at the lower endpoint), 
reach a minimum value of the radial coordinate, and then move up again, approaching 
Release (t) as t — > oo. Notice from (|26|) that, at any given z, the light half-cone i£ ull can 
only point towards increasing z for sufficiently large II (and, for any given II ^ 0, it 
will point upward for sufficiently large z). This means that the entire upward portion 
of ^bh(^) must lie within the region of maximal disturbance of the worldsheet, i.e., in 
the diagonal swath between z grah (t) =t — t gTSih + z m and z re i e ase(£)- 

As usual, determining the exact location of the event horizon is difficult due to 
the global character of its definition: one must know the entire history of the string 
and then integrate (1261) . subject to the stated final condition. It is, however, easy to 
pinpoint with this same equation the location on the worldsheet where z^J^t) = 0, 
which gives a lower bound on the upward portion of zbh(^) (where, by definition, one 
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has i£ u ii(£) > 0). From (j26j) . this happens at 

_ y/l - z 4 m U 2 

2 e rgo(£) — ~ ^ j K^'J 

Z m lL 

where, again, the external force II is meant to be evaluated at the retarded time 
appropriate for the given t, z. At any point along this curve, the downward light half- 
cone is horizontal (or, equivalently, g u = 0), so timelike trajectories must necessarily 
point towards larger z. In other words, z ergo (t) is a stationary- limit curve, and the 
region between it and the upward portion of the event horizon is the analog of an 
ergosphere, a concept whose relevance has been noted previously in the T > context 
in [29]. Setting II = in ([2~?]) implies z ergo — > oo, so, unlike the event horizon, the 
stationary limit curve is located fully within the diagonal region between z gra b(t) and 
z TC \ ease (t). It follows from the definitions of the two curves that z B n(t) crosses z ergo (t) 
precisely when the former attains its minimum value, and so the downward portion 
of the horizon lies below and to the left of z ergo (t). The situation is summarized in 

Fig.m 
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Figure 1: Schematic illustration of the string worldsheet (shaded in gray), in the 
static gauge t = t, a = z, showing the upward null Mikhailov (fixed ttret) lines z gra b 
and Release (solid green), the stationary limit curve 2; er g (dotted blue), and the event 
horizon z B h (thick dotted red) above which lies the worldsheet black hole (shaded 
light red). See text for discussion. 

It is interesting that at T = the notion of a worldsheet black hole plays as much 
of a role as in previous analyses at finite temperature. The appearance of such causal 
structure is seen then to be intrinsically tied to energy dissipation, be it within a 
thermal plasma or in vacuum. We will examine the former case in Section 13.31 It 
would be nice to develop this picture further by exploring the relation between the 
rate at which energy crosses the black hole horizon and the modified Lienard formula 
in (PSJ. 
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3 Single Quark Evolution: Finite temperature 



Having understood the rate of energy and momentum loss and dispersion relation for 
a heavy quark that moves in the SYM vacuum, in this section we restore < oo— 
and consequently h < 1 — in the metric (pQ), to study the same quantities in the 
case where the quark moves through a thermal plasma. In this case, the position 
z = z m < Zh where the D7-branes 'end' is related to the Lagrangian mass m ^> \f\T 
of the quark through [7] 



1 2imi 



z m vA 



8 \ 2m 128 \ 2m / 2m 



{2i 



A thorough generalization of Mikhailov's analytic results [52] to this finite tem- 
perature setup would require finding the exact solution to the Nambu-Goto equation 
of motion for the string on the AdS-Schwarzschild background, for any given trajec- 
tory of the string endpoint at z m > 0. Sadly, we have not been able to accomplish 
this feat. Nevertheless, based on the results discussed in the previous subsection, we 
expect the total energy of the string at any given time to again decompose into a 
surface term that encodes the intrinsic energy of the quark and an integrated local 
term that reflects the energy lost by the quark. 

A priori it might not be obvious that the rate of energy loss in the presence 
of the strongly-coupled non-Abelian plasma should be given by some expression that 
depends just on the behavior of the quark at the given instant, and not on its previous 
history. But in the AdS/CFT context, this property is strongly suggested by the fact 
that the energy of the string is given by a local expression on the worldsheet just 
as much in the AdS-Schwarzschild background that is dual to the thermal plasma 
as in the pure AdS background that corresponds to the SYM vacuum (which, one 
should not forget, is in itself a nonlinear medium). Starting with E written as an 
integral over z, it should again be possible to project back to the boundary along null 
trajectories, to obtain a formula that depends locally on the quark worldline x(t). The 
main difference with the zero-temperature case would be that the null trajectories are 
no longer straight lines. 



3.1 Constant velocity 

We have verified that these expectations are borne out in the case of the only finite- 
temperature solution that is known analytically: the stationary configuration of [7J[B], 



X{t,z) 



t-^ln 



Zh + z 
Zh - z 



Z h , _x 

+ — tan 



(29) 



which describes a quark moving at constant velocity v. Any given point (t, z) on this 
worldsheet is connected to a point (t ret , 0) on the AdS boundary by a null curve t(z), 
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defined by 



(dt\ 4(4 + VT^z 2 ) 
UA ret (4-z*)(z* + zlVT= 



This equation can be integrated to give 

, , Z h (z h + Z \ Z h _! / Z , 



Zh _i 

tan 



^/.(l-v 2 ) 1 /* 

in terms of which the stationary solution 029|) can be written in the form 



tan 



x(t r 



(1-v 2 ) 1 / 4 \z h (l-v 2 )V\ 
Knowledge of this solution allows the total energy of the string, 



E(t) = 
to be reexpressed as 



Vx f Zh , „, Vx r Zh . hx'" + i 



2tt 



dzUl 



2tt 



dz 



z\/l + hX' 2 -% 



Vx 



ret 
00 Z 



2tt 



+ 



Z m \J\ - v 2 z h (l-v 2 )* 

1 



tan 



z h \/\ - v 2 z h (l-v 2 ) 
Similarly, the total momentum of the string, 



z h (l - V 2 )i 



tan 



(1 -v 2 )\ 



Vx 



--h 



Vx 



P(t) = 7— / dzW x = ^— dz 



2ix 



can be rewritten in the form 



2tt 



z 2 hxll + hX' 2 - f 



2n 



2vr 



+ 



z m y/i - 2^(1 - w 2 ) 



tan 



z h (l - V 2 )i 



tan 



z h yjl - v 2 z h (l-v 2 )l l(l-u 2 )4 



As expected, the integrated term in the top line of (1341) and (136|) recovers the 
result for the stationary rate of energy and momentum loss obtained in [TIE], 



The terms in the second and third line of flM|) and fl36l) . then, codify the energy E q 
and momentum p q that are intrinsic to the quark. We can see that, just like at zero 
temperature, dE q /dp q reduces to v only in the pointlike limit z m — > 0. 

An important difference with respect to the T = case analyzed in the previous 
subsection is that here the surface contribution that determines the quark dispersion 
relation arises not only from the lower (z = z m ) but also from the upper (z = z^) 
endpoint of the string. This is in fact the generic situation in the T > case, and 
holds even for the static radial string, where there is of course no energy loss term 
and E is given just by a boundary contribution that defines the thermal rest mass of 
the quark, 



Clearly the second term in (1381) . just like the terms in the third line of (13^|) and 
( I36p . arises from the string endpoint located at the black hole horizon. It is natural 
then to wonder to what extent these terms should be regarded as a contribution to 
the intrinsic energy of the quark, because their value at any given time does not 
depend on the parameters v and F (or d z X(t, z m )) that characterize the state of the 
lower endpoint at the same instant. In fact, since marks the position of an event 
horizon, for any finite coordinate time t the value of the surface contribution at Zh 
is not influenced by the behavior of the z < Zh portion of the string, but depends 
only on the string's configuration at t — > — oo. The same interpretation can be 
then carried over to the gauge theory: the seemingly extraneous terms represent a 
contribution to the energy of the state that depends solely on the initial configuration 
of the quark+plasma system. Throughout the evolution, causality guarantees that 
the behavior of the SYM fields at spatial infinity can only be affected by the initial 
configuration at t — >• — oo, so we can equivalently think of the surface terms at Zh as 
encoding information on the asymptotic boundary conditions for the system. Indeed, 
for dynamical processes, the radial location z = Zh in AdS-Schwarzschild corresponds 
to the deep IR of the gauge theory. 

We conclude then that, to the extent that we wish to compare the energies of 
configurations with different initial/boundary conditions, it is important to keep track 
of the terms arising from the surface contribution at the horizon, despite the fact that 
they are generally independent of the parameters v and F associated with the intrinsic 
dynamics of the quark. Notice, in particular, that with this interpretation the negative 
sign in the second term of (I38p — which is at first sight unexpected when viewed as the 
leading-order thermal correction to the quark mass [T^ — becomes easier to digest: 
it reflects the screening effect of the plasma on the long-range gluonic fields set up by 




(37) 




(38) 
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the quark, which implies a reduction of the energy stored in the IR, in comparison 
with the T = case. 

It is instructive to compare (1341) with the alternative split achieved in [TJ, 

Ed) = ^ ; 2 A * + ^ m _ n 1 , (39 ) 

with Ax = X(t, z m ) — X(t, Zfi)- If one could interpret this latter quantity as the total 
distance traversed by the quark since the beginning of time, then the first term in (1391) 
would give the overall energy lost by the quark, at the known rate ( J37j) . in the total 
elapsed time Ax /v. In view of (1381) . the second term would then imply a standard 
relativistic dispersion relation for the quark, with mass M rest . 

The problem with this interpretation, however (alluded to already in [7J), is that it 
does not properly address the issue of initial conditions. The actual distance travelled 
by the quark is by definition X(t, z m ) — X{— oo, z m ), which agrees with Aa; only if 
X(— oo, Zm) = X(t, Zh)^ This last equality would hold if we had started at t — > — oo 
with the quark at rest (i.e., with the string static and completely vertical), but the 
energy-loss term in (139]) makes no allowance for an initial period of acceleration. 

In contrast with this, the separation obtained in (134"1) has a clear geometric origin 
in the context of the generalization pursued here of Mikhailov's work j52j to the finite- 
temperature case. Within this framework, the portion dE(t, Zh) of the total string 
energy E(t) that is contributed by the segment of the string in the immediate vicinity 
of the horizon (located at X(t, z^) in the notation of (1291) . or X(t ret , Zh) in the notation 
of (1321) ) was lost by the quark at a particular time t ret in the distant past that can be 
deduced from (13TT) . Since this is the highest segment of the string that contributes 
to E(t), t ret marks the precise instant when we need to begin our accounting of the 
energy lost by the quark. At t = t ret the lower endpoint of the string (and, hence, the 
quark) was found a distance d = X(t vet , Zh) — X(t Ict , z m ) behind the location of the 
upper endpoint at t, which according to (1321) translates into 

ZhV _x / 1 \ Z h V _ x 



d = -t— — 7— tan — — -n — — — -n tan 



1-^2)1/4 ^(1 - ^ 2 ) 1/4 / (1-f 2 ) 1 / 4 V^(!-f 2 ) 1/4 

By this logic, the total energy lost by the quark is given by an expression of the same 
form as the first term of (1391) . but with Ax replaced by the actual total distance 
Ax + d = v dt ret . And indeed, we see that the integrated term in the first line 
of (134"1) is larger than the putative energy loss term in (1391) precisely by the amount 
(d/v)dE q /dt, and, correspondingly, the intrinsic energy of the quark identified in the 
second and third line of (1541) is smaller by this same amount than what (13"9"|) would 
have indicated^ 



When writing expressions like X{— oo, z m ), we of course have in mind evaluating the correspond- 
ing quantities at a time that is fixed but arbitrarily far in the past. 

4 Thc fact that, in going from (j34)) to (f39|) . part of the total derivative has been shifted back to the 
integrated term might give the impression that the split between the intrinsic energy of the quark 



20 



3.2 Accelerated quark 

Having gained some intuition from the analysis of a quark moving as in [7J, [8] at 
constant speed relative to the strongly-coupled plasma, let us now turn our attention 
to the more general situation where the quark accelerates. In [TJ, a few important steps 
were taken to have a better understanding of this case; in the next few paragraphs 
we will briefly review the key results. 

A quark that undergoes any type of (forced or unforced) motion (and, in the 
former case, is thereafter released) will be slowed down by its interaction with the 
plasma, and eventually come to rest. The authors of [7J studied the late-time (and 
consequently low-velocity, low-acceleration) behavior of such a quark, by considering 
small, exponentially damped fluctuations around the final rest configuration. In dual 
language, this involves a determination of the quasi-normal modes on the worldsheet 
of a static and purely radial string. From their analysis they were able to numerically 
deduce, for any given quark mass parameter z m , the value of the drag coefficient 

A* - --^ • (40) 

P Pq dt V ' 

As long as one maintains the restriction to the non-relativistic regime, the above 
definition is equivalent to \i — —(l/v)dv/dt, and by construction yields a result that 
is independent of p (or v). A few representative values were tabulated in [7J. 

Additionally, the authors of [7J gave an analytic derivation of the low-velocity 
dispersion relation for the quark, which they found to take the form 

£ ' = M '» + 2|- + 0( *' (41) 
where M rest is the thermal rest mass fl38|) . and 

A4,„ - ^ (42) 

the kinetic mass of the quark. In the heavy quark limit m ^> a/AT (z m -C z/j), where 
([38D and (TJHD imply that 

M M = m -^ + oU^ | , (43) 
they found from their quasi-normal mode calculation that 



fi = 



2 m 



2m 




(44) 



and the energy that has already been lost is inherently ambiguous. The wide latitude available in 
this case, however, stems from the steady-state nature of the configuration under scrutiny. In the 
general case, clearly it is a very non-trivial property for a particular contribution to the energy E{t) 
of the string to be expressible as a functional only of the state of the endpoints at the given instant. 
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which through (142]) leads to 



M kin = M rest + O 



in 



2m 



(45) 



It was noticed in [7j that expressions ( 14TT) and (1441) . valid in the low- velocity, low- 
acceleration regime, as well as the value of the drag coefficient ( 1401) deduced from 
( |37|) . valid for a quark with constant but otherwise arbitrary velocity, are consistent 
with a relativistic dispersion relation of the form 



E q = M rcst - M kin + J Pi 



ML 



M rest + M kin ( 7 - r 



(46) 



Additional evidence for this relation was found from the analysis of back-to-back 
motion of a quark and antiquark formed within the plasma. The results of [7J in this 
setting will be reviewed and extended in Section HI 

In view of the discrepancy between 



61) and the relation E„ 



that 



one would naively infer from the second term in fl39l) . the authors of [7J emphasized 
the need for a more detailed study of the quark's intrinsic dynamics, and proposed 
a plan of attack. They observed that if one starts with the quark at rest in the 
hot medium, and then accelerates it with an external force, then as long as energy 
dissipation is negligible, it is natural to define the total intrinsic energy of the quark 
as the initial rest energy plus the work done by the external agent. Motivated by this 
proposal, we have studied the early-time behavior of a quark initially at rest, which is 
accelerated by a time-dependent external force that is turned off after a short period 
of time, allowing the quark to move thereafter only under the influence of the plasma. 

As we already mentioned, regrettably, we have not been able to find a general exact 
solution to the Nambu-Goto equation of motion for the string in the non-stationary 
case. Nevertheless, it is certainly possible to find numerical solutions that describe 
string configurations dual to the gauge theory setup described above. Starting at 
t — 0, a string extending from a fixed value of the radial coordinate z = z m to the 
horizon at z = Zh, initially at rest and vertical, is accelerated by applying an external 
force F(t) to its lower endpoint. After a time ^release, this external force is set to zero 
and the string moves freely in the curved background. 

Using (pQ) and ([2]) and imposing the condition that the string moves only in the 
x = x 1 direction, the equation of motion is 



d_ 

dz 



-hX' 



1 + hX ' 2 - ^ 
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z 2 h\ I 



hx' 2 - 4^ 







(47) 



and we must set the initial and boundary conditions to be 

X(0,z) = 0, X(0,z) = 0, X'(t,z m ) = fit) , X(t,z h ) = 0. 



(At 
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The first two conditions here simply implement the requirement that the string start 
out being static and purely radial. In the third condition, for a given external force 
F(t) = (vA/27r)II^(£, z m ) acting on the quark, the function f(t) could be determined 
using the relation ([3]) between X' and 11*. In practice we find it easier, however, 
to specify f(t) and use the results of the numerical integration together with ([3]) 
to deduce the associated F(t). By studying a number of different examples, we 
have verified that the instantaneous rate of energy loss after the quark is released is 
independent of our choice of f(t) and £ re i e ase, i-e., the quark does not care about its past 
history For the trajectories that we will plot below, we have used f(t) = bt(t—t TC i ease ), 
with t re iease = 0.3 / tcT and adjustable b. 

The fourth and final condition in (1481) specifies that the string endpoint at the 
horizon remain fixed, reflecting the fact that the wavefront for the disturbance pro- 
duced by the external agent at the z = z m endpoint of the string will not reach z = Zh 
before an infinite amount of (boundary) time has elapsed. In order to implement this 
condition in our numerical integration, we have set X(t, z maiX ) = at a radial cutoff 
z-max — 0.999^, and considered in all cases an integration time t max smaller than the 
time /_ s f max dz/(l — {zj ' Zh) A ) that it takes the wavefront to reach the cutoff. 

For applications of this formalism to phenomenology, we must choose values of 
the mass parameter z m based on the charm and bottom quark masses, m ~ 1.4,4.8 
GeV. The issue of how best to translate between the SYM and QCD parameters has 
been discussed in [48]. Taking oiqcd = 0.5 (gQCD = v^vr), N c = 3 and Tq C d = 250 
MeV, and employing the "obvious" prescription gyu = 9qcd and Tsym = Tqcd, 
from (1281) we find that z m /zh ~ 0.40 for charm and z m /zh ~ 0.11 for bottom. If, on 
the other hand, one uses the "alternative" scheme gy M iV ~ 5.5 (motivated in [48] 
through a rough matching of the AdS/CFT and lattice quark- ant iquark potentials) 
and Tsym — 3~ 1//4 Tqc£> — 190 MeV (which follows from equating the energy densities 
of the two theories), then (l28~l) leads to z m /zh ~ 0.16 for charm and z m /zh ~ 0.046 
for bottom. In our analysis, we have covered a significant range of masses, but below 
we will present the results only for three representative values in the neighborhood 
of the charm mass: z m /zh = 0.2, 0.3, 0.4. 

We have carried out the numerical integration of (1471) subject to (148]) using the 
NDSolve routine of Mathematica 5.2. Based on the variation of our results upon 
doubling the number of integration steps, we estimate our numerics to be accurate 
to better than 1%. The integration time shown in all plots is given in units of 1/ttT, 
which (for Tqcd = 250 MeV) corresponds to 0.25 fm/c under the "obvious" and 
0.33 fm/c under the "alternative" prescription of [38]. Unfortunately, the numerical 
integration degrades rather quickly, so in either scheme our investigation is limited 
to intervals that are an order of magnitude below the experimental timescale of the 
plasma (typically ^breakdown ~ 0.9/7rT ~ 0.3 fm/c). For the same reason, even though 
the quark can be taken to relativistic velocities at the point of maximal acceleration, 
we can only achieve rather small velocities at the time of release (for the most part, 

^release *^ 0.1). 

Our results show a qualitative difference between the initial stage (0 < t < t Te i ease ) 
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where the quark is accelerated by means of the external force F(t), and the second 
stage (t re i e asc < t < t max ) where it moves only under the influence of the plasma. We 
will begin by discussing the latter stage, which would appear to be more relevant 
from the phenomenological perspective. 




Figure 2: Quark velocity as a function of time from our numerical integration (in red) 
compared against (J5T]) with the value of \x deduced in [7J (in black), and floTT) with fi 
chosen to fit the data (in light blue), for a) z m j z^ = 0.2 and b) z m /zh = 0.4. See text 
for discussion, 

The most direct way to inquire whether the output of our numerical integration 
for the accelerated quark conforms to the the constant-velocity (or late-time) results 
of E] is to compare the corresponding quark trajectories. If one assumes the 
dispersion relation (I46p . then the equation of motion for a quark subject only to the 
drag force fj4*0|) with constant \i is 

| = (49) 

whose solution is [7J 

v(t) = Vvdcasc = . (50) 

) 2 . -I- ( I — ?; 2 , ")f>2/i(t— Release) 
release 1 V release/ 

As we already mentioned, in our numerical results ^release € 1, so we are only able to 
test the non-relativistic version of (I5DD . 

V(t) = Wsee-^-'-l-se ) (51) 

(and, since our integration is limited to small time intervals, we would in effect see just 
the linear portion of this function). A comparison between this analytic prediction 
and our numerical results for t > t TC \ ease is given in Fig. [2j It is evident from this plot 
that, in the early stage of motion covered by our analysis, the quark dissipates energy 
at a rate much lower than the late-time result of [7J. Indeed, for z m / z h = 0.2, 0.3, 0.4 
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the asymptotic friction coefficient is respectively /xi ate /7rT = 0.25,0.41,0.59, but our 
numeric results for v(t) are best approximated by /icariy/ 71 "^ = 0.08,0.15,0.26. 

We can also attempt to perform the comparison directly at the level of energy loss 
rates. An important drawback of working with the numerical solution, however, is 
that we cannot achieve a direct splitting of the total energy E of the string, as we did 
in the stationary case (as well as in the general case at zero temperature). This means 
that, a priori, we know neither the correct form of the quark dispersion relation nor 
the formula for the rate of energy loss. But, given that their sum remains constant 
throughout the evolution, finding a prescription for one of these two quantities would 
enable us to compute the other. 

The authors of [7] assumed that the energy loss would be negligible for sufficiently 
short acceleration intervals (i.e., for small t re iease), because the drag force exerted by 
the plasma would not have been able to perform a substantial amount of work. If 
true, this would allow a direct empirical determination of the dispersion relation. 
Unfortunately, the situation is not so simple, because, as we learned in Section [21 
the quark loses energy through radiation even in the absence of the plasma, and this 
effect must be taken into account to establish what fraction of the total string energy 
E is intrinsically ascribable to the quark. 

To attempt to cut this Gordian knot, we should recall, from our study of the cases 
where we had analytic control, that the dispersion relation arises as a surface term, 
with contributions from both endpoints of the string. Based on our previous results, 
we expect the dispersion relation for the quark in the thermal medium to take the 
form 

J' Z ^ 4 tt2 T - -) + • ( 52 ) 

y/(l-V 2 )(l-ziU 2 ) ZhJ 

The first term here has been copied from the T = expression (1191) (again abbreviat- 
ing I! = II*), and is meant to approximate the contribution from the string endpoint 
lying on the D7-branes. The second term arises from the endpoint that reaches the 
horizon, which we understood above to encode information about the initial con- 
ditions. Knowing that our starting configuration is static, we can simply read off 
this contribution from the second term in the rest energy (]38p . which we have in- 
terpreted as a screening effect. Including this term is important to ensure that (|52[) 
reduces to the correct result in the case of a quark that is static and unaccelerated, 
E q (0,0,T) = M rcst . 

Clearly equation (152]) is just an approximation, because the first term should 
receive thermal corrections. In particular, it is natural to expect the factor of 1 — z^Jl 2 
in the denominator to be replaced by h(z m ) — z^JI 2 , so that, just like in the T = 
case, the quark energy diverges at the critical value of the electric field, which now 
corresponds to II = h(z m ) / z^ [20]. More generally, since the background metric 
(JTJ) knows that T > only through the factor h < 1, any corrections to the z = z m 
surface term due to the presence of the plasma should be of order 1 — h or h', and are 
consequently small in the heavy quark regime z m j z/,< 1. For the phenomenologically 
interesting values z m /zh = 0.2,0.3,0.4, the error in the dispersion relation ff52]) is 
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estimated to be ~ 4 — 16% (much larger than our ~ 1% numerical error). 

Notice that (1521) without any corrections should give the exact dispersion rela- 
tion in the infinite-mass limit z m — > that has been the focus of many AdS/CFT 
investigations of energy loss (e.g., [8] [25]). In this limit one is left with 



i.e., the F-dependence drops out and, just like in the T = context at the end of 
Section[2, one recovers pointlike behavior. We can see that, as expected, the first term 
in (|53|) agrees with the z m — > limit of the second line of the stationary expression 
(134]) . but the second term in (|53|) disagrees with the limit of the third line of (134j) . 
due to the different initial conditions. 

Expression (152"]) as a whole looks superficially rather different from the dispersion 
relation (T4"6"|) proposed in [7]. In particular, (T4"6"|) evidently cannot reproduce the F- 
(or II-) dependence seen in fl52l) . whose presence is supported by the zero-temperature 
results of the previous subsection. It is important to remember, however, that such 
dependence indeed would not have been visible in the quasi-normal mode analysis 
used to derive (14T1) or in the quark-antiquark evolution that gave part of the support 
for ( j46|) . because in those calculations the external force was taken to vanish. 

Setting n = 0, ([52]) reduces to E q = ( v / A/2tt)(7 / z m - l/z h ) + 0(z 2 m /zl), while 
f ]46]) translates (via (EHJ) and (ED)) into E q = (V\/2tt) {^/z m - y/z h ) + 0(z m /z 2 h ). 
The two expressions differ in the form of the second, z m -independent term, which we 
have understood to encode the initial conditions for the gauge system. Given that the 
initial conditions for all situations considered in |7] differ from our current setup, there 
is no reason why we should expect the corresponding terms to agree. Beyond this, 
there also appears to be a discrepancy in the order of magnitude for the corrections 
to the two expressions. We remain puzzled by this apparent mismatch, because as 
explained below (I512"]) . we do not see how the zero-temperature relation could receive 
corrections higher than order 1 — h or h'. 

We interpret the w-dependence seen in the surface contribution at Zh implied by 
( 146]) to be a reflection of the fact that, for the quasi-normal mode considered 
the string endpoint at the horizon is indeed moving, unlike what happens in our 
In other words, we find that the 'kinetic mass of the quark,' defined as the coefficient 
of the v 2 /2 term in E q , is sensitive to the initial conditions. In particular, for a quark 
that is initially static, we have obtained a dispersion relation of the same relativistic 
form as fj4"6"]) . 

E q = M rest + M kin (7 -l) + \ m( j | , (54) 




5 For the quark-antiquark configurations that were analyzed in [7] and will be reviewed and 
generalized in Section 31 there is no endpoint at the horizon, so we should expect to get only a 
surface contribution from z = z m . Independently of that, the evolution in that case is not sensitive 
to the value of M\^ n , which appears on both sides the equation of motion (|75|). and consequently 
drops out, leading to ([49]) . 
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but with M kin = \ r \/2irz m = m + 0(m(y/\T/2m) i ). 

We are finally in position to address energy loss. For the type of process we 
consider, the total energy of the string at time t is given by 

E(i) = M rost + E input (t) , 

where the second term is the work performed by the external force, which can be 
computed by integrating the energy current — Ilf (t, z m ) shown in ([3]). The energy lost 
by the heavy quark can then be obtained as 

E lost {t) = E(t) - E q (v(t), E(t),T) = E mput (t) - E kin {t) , (55) 

where the last term is the kinetic energy 

£kin = EJv, F, T) - M rcst = ——( - 1 ~ Z ™ vU - 1 ) . (56) 

Of course, in the t > t re \ ease (II = 0) stage, with the small velocities that we achieve 
we are only sensitive to the non-relativistic (quadratic in v) terms in this equation. 

E E 
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Figure 3: a) Accumulated energy loss (in units of yXT/2) as a function of time (in 
units of 1/vrT) for a) z m = 0.2 (red) and z m = 0.4 (blue) with f rc i C asc = 0.051, and b) 
z rn = 0.2 (red) and z m = 0.3 (blue) with f re iease = 0.073. For comparison, the dashed 
curves of the same colors give the energy loss that would follow from the stationary 
rate (1571) obtained in [7J [S]. The green lines represent the rate (1571) evaluated with 
v = v release , which can be contrasted against the slope of the numerical curves, shown 
in black. 

In Fig. [5J we compare the energy loss at for different masses but equal velocity 
at the time Release when the external force is set to zero. The curves have essentially 
constant slope for t > t re \ ease , meaning that the rate of energy loss is nearly constant 
in the limited time window that we have access to. The lines in the figure contrast 
the instantaneous rate of energy loss in this interval, obtained by a numerical fit and 
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denoted henceforward by (d t E q ) n , against the corresponding stationary result of [3 [8], 
which we will denote by (d t E q ) s . For z m = 0.4, (d t E q ) n is close to two times bigger 
than (d t E q ) s , and as the value of z m decreases (corresponding to heavier quarks), the 



difference between the two rates grows bigger. In particular, for z„ 
at least three times bigger than (d t E q ) n . 



0.2, (d t E q ) s is 
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Figure 4: a) Accumulated energy loss (in units of \f\T '/2) as a function of time (in 
units of 1/vrT) for a) z m = 0.2 with v = 0.056 (red) and v = 0.111 (blue) and b) 
z m = 0.3 with v = 0.036 (red) and v = 0.073 (blue). The dashed curves of the same 
colors give the energy loss obtained with the stationary rate (I37j) . The green lines 
represent this rate evaluated at t = t Te \ ea , se with velocity v, which is to be contrasted 
against the slope of the numerical curves, shown in black. 

A similar comparison is shown is Fig. HJ but with a fixed mass value and different 
release velocities. In all cases, we have found again that the rate of energy loss (]37p , 
valid in the stationary regime, is above our numerical result. The discrepancy ob- 
served is significantly larger than our estimated margin of error. Our results therefore 
provide clear evidence that, in the mass range of primary phenomeno logical interest, 
there exist conditions under which the rate of energy loss for a heavy quark that 
moves only under the influence of the plasma can be substantially smaller than the 
rate obtained in [7] for the steady-state or late-time configuration. We will return to 
this point in Section 15.41 

As one would expect, the numerical rate of energy loss depends on z m and v. As 
shown in Fig. [SJ for fixed mass values, it varies quadratically with the velocity of the 
quark. As we had mentioned before, for the small velocities that we can attain, we 
must deal only with the non-relativistic approximation of the dispersion relation (151)1) . 
This implies that the accumulated energy loss (155]) varies quadratically with the quark 
velocity v, and so the constant slope that we are reading off of the t re iease < t < t max 
portion of the numerical curves in Figs. [3] and H] should be proportional to vAv, where 
Av denotes the small change in velocity in the given time interval. The parabolic 
behavior seen in Fig. [51 then, tells us that Av oc v. In other words, in the leading 
non-relativistic version of (T49l . the friction coefficient /i is independent of the quark 
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velocity, as one would expect. 

Unfortunately, due to limitations with the numerical integration, we have not 
been able to characterize the dependence of the dissipation rate on the mass of the 
quark, beyond the statement that (dtE q ) n increases roughly linearly with increasing 
z m /zh (decreasing m). In particular, above z m /zh ~ 0.75, the numeric rate becomes 
indistinguishable (within our margin of error) from the stationary result fl37|) . 



(d t E q ) n (d t E q ) 
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Figure 5: Rate of energy loss (points, in units of a/A7tT 2 /2) as a function of velocity 
for a) z m = 0.2 and b) z m = 0.3, together with the corresponding quadratic fits 
(d t E q ) n (0.2,v) = 0.31v 2 and {d t E q ) n {0.3,v) = 0.4k; 2 . 
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Figure 6: Comparison between the accumulated energy loss (in units of y/XT/2) 
versus time (in units of 1/vrT) using ( |55i) in blue and the modified Lienard formula 
( 1T81) in red for a) z m = 0.2 and b) z m = 0.3. See text for discussion. 

Let us now consider the initial stage < t < t re \ ease . From the beginning portion 
of the curves in Figs. [3] and HJ we can see that the situation when the quark is 
subjected to an external force is opposite to what we described above for unforced 
motion: the rate at which energy is dissipated can be substantially larger than ( |37l) . 
suggesting that in this segment of the quark trajectory the mechanism of energy loss 
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is qualitatively different. In particular, the energy lost in this region is by no means 
negligible, which is precisely the obstacle that prevents us from directly inferring the 
quark's dispersion relation from our numerical results, as envisioned in [7J. We had 
anticipated this already in the discussion above, based on our results for the T = 
case in Section [21 

When, using the numerical data of the evolution within the plasma, we compare 
as in Fig.[H]the energy lost by the heavy quark, Eq. ([51)1) . against the modified Lienard 
formula ([TBI) , we find that for t < t ie \ ease the two curves are virtually indistinguishable 
from one another when z m /zh = 0.2, 0.3, and, to a lesser extent, 0.4. In fact, this 
approximate agreement continues to hold (albeit somewhat reduced) even if in (fl8l) 
we plug in the data of the quark trajectory at zero temperature. In other words, for 
quark masses in the phenomenologically interesting range, the quark behaves initially 
as if there were no plasma, and loses energy through radiation, at a rate equal to the 
modified Lienard formula given by the first term in fTTHl) . This is of course as one 
would expect from the gauge theory perspective, for a heavy quark should indeed be 
insensitive to the plasma for very early times. 

On the gravity side, the issue is that, for these relatively low values of z m /zh, the 
factor of h is so close to unity that the modifications induced by the black hole horizon 
on E inpnt and E\ ost only become appreciable when a sufficiently large time has elapsed. 
We should emphasize that it is not possible to reproduce the numerical results for 
v (t) and E input (t) = E^ in {t) + E\ ost (t) in this initial stage using the dispersion relation 
( )46l) and dissipation rate (j37l) derived in [7J (the values of the latter are indicated by 
the dotted curves in Figs. [HH). 

Beyond t re iease the curves in Fig. separate, showing the influence of the hot 
medium, although somewhat diminished compared to the stationary result fl37|) . as 
we know from Figs. [3HJ Since it was shown in [7J that (137|) will hold at asymptotically 
late times, we expect the rate of energy loss seen in Figs. [2][6]to increase as the system 
evolves further. Regrettably, with our very limited integration time we are not able 
to track the evolution far enough to locate the characteristic transition time to the 
asymptotic behavior. We will return to this issue from a different perspective in 
Section WM 

3.3 Late-time behavior and worldsheet black hole 

It is interesting to visualize the evolution of the system beyond the limited time 
interval covered by our numerical data, in parallel with our discussion for the zero- 
temperature case in Section 12.31 The quark, initially static, and accelerated by an 
external force F(t) between t = t grab (originally t = 0) and t = t re \ easc , will thereafter 
decelerate under the influence of the plasma, approaching rest at some location as 
t — > oo. In the dual gravity description, this means that the final string embedding, 
just like the initial (t < t gra b) one, must include a static vertical segment extending 
all the way from the D7-branes at z = z m to the black hole horizon at z = Zh, to 
represent the quark at rest. 
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As time progresses, the lower (z = z m ) string endpoint traces out the trajectory 
of the quark, moving from x = to x = The upper (z = Zh) endpoint, on the 
other hand, remains at x = for all finite times, because the wavefront generated 
on the string by the acceleration of the bottom tip will reach the spacetime horizon 
only at t — > oo. Given these boundary conditions, on the x-z plane the string will 
clearly evolve from purely vertical to -i-shaped, with a horizontal segment at z — z^, 
extending from x = to x = x^. All of the energy (and momentum) lost by the quark 
throughout its evolution ends up in this top portion of the string, which encodes the 
IR region of the gauge theory The region of no escape is bounded as in Section 12.31 
by a worldsheet horizon that can be identified by standing on the (£oo; Zh) corner of 
the string and projecting back along the downward-pointing light half-cone z^J { (t), 
now defined by 



The resulting curve Zbh(^) descends from z — > oo before t = t gra b; reaches a minimum 
value of the radial coordinate, and then moves up again, finally approaching the 
'corner' as t — > oo. The upward segment lies fully within the maximally-disturbed 
region of the worldsheet, located between the upward null curves z gTa b(t) and z re \ ease (t) 
obtained by integrating (l57j) with the upper choice of sign and with initial condition 
"2 g rab(£grab) = z m or Release (Release) = z m , respectively. A lower bound on the location 
of the upward portion of zbh(^) is given by the (upward segment of the) stationary 
limit curve z ergo (t), defined as the locus where g tt = 0, i.e., X 2 = h. The situation is 
summarized in Fig. [7J 

If, instead of releasing the quark, we pull it with constant velocity v for an arbi- 
trarily long period of time (i.e., if t re i ea sc — > oo), then we approach the steady-state 
configuration of [5] , and in so doing stabilize the worldsheet horizon (as well as the 
stationary-limit curve) at z = z v = (1 — v 2 ) 1 ^, just like in [2U1 US] (and [29J ) . 

The behavior of the string long after t re iease was determined quantitatively in [7J, 
and is of the form X(z, t) = Xoo — A(z)e~' J ' latet , with the friction coefficient //i a t e given 
by the lowest quasi-normal frequency. It was shown there that the imposition of a 
purely ingoing boundary condition near the spacetime horizon z = Zh forces the string 
to deviate from the vertical by a divergent amount (A(z) — > oo as z — > z^). This is 
consistent with the development of a 'corner', just as we have argued above for the 
case where the quark is initially static. 

Just as we found at T = 0, we see here that energy dissipation seems to be 
irrevocably tied to the appearance of a worldsheet black hole. Since our numerical 
results show that the initial energy loss is controlled by a friction coefficient /! ca riy < 
Aviate; h would appear like in the unforced case the asymptotic rate sets in only when 
the string is sufficiently close to its final -i shape. It would be interesting to establish 
a more detailed connection between the instantaneous rate of energy dissipation from 
the quark and the rate at which energy crosses the worldsheet horizon, generalizing 
the results of [2"Ul [TU] for the stationary configuration. 




(57) 
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Figure 7: Schematic illustration of the string worldsheet (the rectangle shaded in 
gray), projected for convenience onto the spacetime x-z plane. To aid the visualization 
of the evolution, snapshots of the string are given for three different instants: at (any 
time up to) t = t gTab (thin dotted orange), when the string is at rest and vertical at 
x = 0; at t = Please (thick solid orange), when it has already been partially deformed 
by the application of the external force F(t); and at t — > oo (thin solid orange), when 
it has adopted a -> shape, and its vertical segment has come to rest at x — Xqo. The 
diagram additionally shows the upward null (fixed t t rct) curves z gra b and Release (solid 
green), the stationary limit curve z eTgo (dotted blue), and the event horizon zbh (thick 
dotted red) above which lies the worldsheet black hole (shaded light red). See text 
for discussion. 
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4 Quark- Antiquark Evolution 



A string with both of its endpoints on the D7-branes describes a quark-antiquark pair. 
The situation that is closest to modeling the dual process of primary phenomenological 
relevance, where a heavy quark and antiquark are created within the plasma at time 
t — and then separate from one another, is such that the string endpoints start out 
with coincident positions but different velocities. In a first exploration of this system, 
it is interesting enough to consider the simple case where the string endpoints are 
taken to separate back-to-back with the same initial speed t>o, meaning that the pair's 
center of mass frame coincides with the plasma rest frame. We will let x denote the 
direction of motion. 



4.1 Review of earlier results 

A numerical study of this problem was carried out in [7] . It was found that the quark 
and antiquark trajectories can be more efficiently followed to later times if instead of 
describing the string embedding X M (r, a) in the obvious static gauge r = t, a = z, 
one astutely chooses worldsheet coordinates for which the constant r slices manage 
to reach larger values of x near z = z m while staying away from the horizon at 
z = Zh- This is most easily implemented by working with the Polyakov (rather than 
Nambu-Goto) action 



S, 



1 



Ana' 



dr / da^g ab G, u d a X»d b X» 



R 2 
2ira' 



d 2 a Cp, 



(58 



with g a b the intrinsic metric on the string worldsheet, and making the non-standard 
gauge choice g TT = —s, g aa = 1/s, g UT = 0. The 'stretching factor' s = s(cr, r) 
implicitly defines our choice of worldsheet coordinates, and is meant to be adjusted 
by hand to keep the numerical integration away from the horizon for the longest 
possible time. For a given trajectory, it is the function z(t, a) that determines the 
proximity to the horizon, so this goal is achieved by setting s = s(z(t, a)) (at the level 
of the equations of motion). Following [7] we will use s oc (zh — z) p in the examples 
below. 

In this setting, the evolution of the string is controlled by the equations of motion 
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(where ' = d T , ' = d a ,) supplemented with the constraints 

- hit' + xx' + h~ l zz' = , (60) 
- h{ ? + s H l2 ) + {x 2 + s 2 x l2 ) + h- l {z 2 + s 2 z' 2 ) = 0, 

which as always amount to the statement that the intrinsic metric on the worldsheet 
must be proportional to the induced metric (thereby establishing the classical equiv- 
alence with the Nambu-Goto formalism). Given initial data that satisfy (160]) . the 
requirement that the constraints continue to hold throughout the evolution gives an 
important consistency check on the numerical integration. The results we will report 
throughout this section were obtained with Mathematica 5.2's NDSolve integration 
routine, with the constraints typically satisfied at the 10~ 5 or 10~ 6 level. 

In terms of the momentum densities IT^ = dCp/d(d a X^ L ), we see that, as usual, 
the first two equations in (159]) express the conservation of the Noether currents IT" 
and IT", respectively associated with invariance under translations in t and x. To 
describe a quark-antiquark pair that is not acted upon by any agent other than the 
plasma, we must choose the standard Neumann/Dirichlet boundary conditions 

t'(r, 0) = t'(j~, 7r) = , x'(t,0) = x'(t,tt) = , z(r, 0) = z(r, tt) = z m Vr . (61) 

The total energy 



E = ^ Tda {-W t ) = — r da 4 (62) 



and x-momentum 

r 2 r , m v 7 ! r r h± 



of the string are then conserved. 

For the problem at hand, the authors of [7] identified a one-parameter family of 
initial conditions that correctly satisfy the constraints ( !60l and are compatible with 
the boundary conditions (|6"T|) . Working from now on in units where z^ = 1/txT = 1, 
these conditions take the form 

(64) 



t(0,a) = , 


i(0,a) 


= A, 


x(0,a) = , 


x(0,a) 


= A\/ 1 — zf n cos cr 


z(0, a) = z m , 


i(0,a) 


= A[l — z^} sin a , 



and describe a string that is pointlike at t = and grows for t > as a result of its 
non-zero initial velocity 



x 



v x{°) = j = V 1 - z^cosa , v z (a) = j = [l-zjsina . (65) 

The parameter A (which has been rescaled here by a factor of a/1 — ^ with respect 
to [7]) controls the energy (1621) of the configuration, 

\/A f, / hi \ y/X(l- z A m )A 



E = / da = -t P \ T \ ■ (66) 
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As expected, the total x-momentum vanishes. 

Using the initial conditions it was found in [7J that, depending on the value 
of E, the subsequent behavior of the string endpoints can be of two different types. 
When the energy of the pair is large enough (essentially, E > 2M rest ) the quark and 
antiquark move apart and are able to escape from one another's influence, so they 
simply slow down monotonically until they are finally (at t = oo) brought to rest by 
the plasma. For low E, on the other hand, the mutual attraction of the quark and 
antiquark manages to stop them and make them reverse direction, after which they 
undergo a number (larger than one quarter) of oscillations before dissipating all of 
their energy to the plasma. During consecutive half-cycles the body of the string is 
alternately above and below the z = z m line, as a result of which the corresponding 
motion is asymmetric J§ 

4.2 Generalized initial conditions 

Intuitively, it should be possible to generalize the initial velocity profiles (1651) pro- 
posed in [7J to more general functions v x (a), v z (a), which amounts to stipulating 
that x(0,a) = Av x (a), i(0,a) = Av z (a). Based on the symmetry of our problem, 
for simplicity we restrict attention to functions v x (a) that are odd on the interval 
[0,7r]. Compatibility with the Neumann boundary condition in (I6TI) requires that 
v' x (0) = 0, the Hamiltonian constraint in (160]) determines v z (a) in terms of v x (a), and 
the Dirichlet boundary condition in (IBTj) demands that v z (0) = 0. Altogether, then, 
we find that we are allowed to choose 

t(0,a) = 0, i(0,a)=A, (67) 
x(0, a) = , i(0,a) = Av x (a) , 
z(0, a) = z m , z(0, a) = Ay/l - z^y/l -z^- v x (a) 2 , 

with v x (0) = a/1 — z^ and v' x (0) = 0, which constitutes an infinite-parameter gener- 
alization of the initial conditions ( 1641) . 

The energy of all of these configurations is given by the same formula ( |66l . Since 
the string is initially a point, it might seem peculiar that one can physically distinguish 
among various possible velocities for its different 'internal points' — one might suspect 
that the seemingly different initial profiles (1671) are all related to each other by gauge 
transformations. That this is not the case can be seen most straightforwardly by 
evolving a few of these initial configurations forward in time: as shown in Fig. [SJ the 

6 Configurations where the string lies below its endpoints have been considered previously in [29l 
130] . In such circumstance, one suspects that the string would prefer to shrink by sliding its endpoints 
along the D7-branes, moving them closer to the boundary and to each other. To examine this 
question, however, one must remember: first, that the boundary conditions for the string are purely 
Neumann or Dirichlet only when the D7-brane embedding is described in the original Cartesian 
coordinates, and are actually mixed in the spherical coordinates that are naturally employed after 
taking the AdS/CFT limit; second, that the string coordinates on the S 5 describe the internal SU(A) 
degrees of freedom of the quark, which one may or may not wish to fix externally. 
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resulting spacetime trajectories are found to be distinct. Notice in particular that 
in the dashdotted trajectory the quark actually reverses direction before coming to 
rest, which shows that the boundary between oscillatory and non-oscillatory behavior 
depends strongly on the way in which the gluonic field (dual to the string) is excited. 
So the moral of the story is that, while the string is initially a point in spacetime, it is 
most definitely not a point in phase space, and it is this fact that allows the existence 
of a truly infinite-dimensional family of initial conditions. 

v x (a) x 



Figure 8: (a) Three different initial velocity profiles for the string: cos a (solid), cos 3 a 
(dashed), and cos3a (dashdotted). (b) The corresponding trajectories for the a = 
string endpoint, for quark mass parameter z m = 0.2. Even though all 3 configurations 
have the same energy E/2M Test = 2.45 and initial quark velocity Vq = 0.9992, the 
evolution of the quark is clearly rather different in each case. 

The question of the gauge-dependence of our description does however serve to 
highlight a useful point: instead of sampling different initial conditions by varying 
the functional form of v x (a) in (1671) (for a fixed choice of the worldsheet coordinate 
a), we can keep the form of v x (o~) fixed and change the meaning of a. One way 
to do the latter is to modify our choice of the initial stretching factor s(0,cr). It 
is easy to see that this indeed leads to different spacetime trajectories, even if the 
initial energies (|6"6"|) are appropriately matched. Conversely, when tweaking s(r, cr) to 
find the choice that optimizes the numerical integration for a given physical string 
configuration, one must keep the initial stretching factor s(0, a) fixed, to ensure that 
the initial conditions for the evolution do not change. For our calculations we found 
it convenient to use s(r, cr) = (1 — z(r, a)) p /(l — z m ) p_1 , with adjustable p. 

It is natural to wonder what the interpretation of the different initial conditions 
(|67|) is in the dual SYM language. The answer is provided to us by the standard 
recipe for correlation functions [2]: different time- dependent string profiles corre- 
spond to different time-dependent configurations of the gluonic fields [56j SI] • Just 
like the string embedding is not completely characterized by giving the location of 
its endpoints, the initial state of the gauge theory is not uniquely characterized by 
specifying the position of the quark and antiquark. This is of course true already 
in zero-temperature QED, but in that case the linear character of the equations of 
motion makes it easy to identify, for a given field configuration, the portion that is 
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directly ascribable to the sources of interest. 

Notice from (I6T1) that the initial velocity of the quark, v = v x (0) = a/1 — z^, 
is not a free parameter of the system, but is uniquely fixed by the choice of z m , or 
equivalently, by the quark's Lagrangian mass m, according to fl28|) . The reason for 
this is easy to understand on the string theory side. At t — 0, the pointlike string 
happens to obey not only i = 0, in compliance with flBTj) . but also z' = 0, according 
to flBTl) . as would befit an endpoint that is free in all spacetime directions. It is 
well-known that the endpoints of such a string must move at the speed of light (see, 
e.g., [61J), and, given that the endpoints are located at z = z m , we see from fll)) that, 
indeed, a coordinate velocity v x = a/I — z^ corresponds precisely to a proper velocity 



V x = v x / y/—G a = 1, a fact that was first pointed out in [35] Jj In the gauge theory, 
this identification confers then a special status to the mass-dependent velocity 



- 4 , (68) 

whose meaning will be discussed further in the next subsection. 

At least from the gauge theory perspective, one would expect to be able to find 
configurations in which the initial quark velocity vq is freely adjustable. Given the 
discussion of the previous paragraph, we see that on the AdS side this can be achieved 
with initially coincident string endpoints only if we choose z' ^ 0. To satisfy the first 
constraint in (160]) we must then set z = 0. Picking for the string a velocity profile 
v x (a) that is an arbitrary odd function on the interval [0, ir], the complete second set 
of allowed initial conditions is then 

t(0,(r) = 0, i(0,a)=A, (69) 
x(0, a) — , x(0, a) = Av x (a) , 

*(0,<7)=C(<7), Z(0,<7)=0, 

with C( a ) an even function on the [0, tt] interval, which satisfies £(0) = z m for com- 
patibility with the Dirichlet boundary condition in (16TI) . and 



to comply with the Hamiltonian constraint in (1601) . By construction, the initial quark 
velocity can now be chosen arbitrarily, as long as Vq < v m in order for the right-hand 
side of (1701) to be real at o = 0, ir. The energy (16"2|) of the configuration is now given 
by 

2t Jo V»*V_o 2t Jo s (0C K 



7 The same reasoning in fact applies to 'all points' on the string: their proper initial velocities 



V *( a ) = n ' V ^ a > = i IT = 



clearly satisfy V% = 1, independently of the choice of v x (a). 
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The linear x momentum is still zero. 

Conditions (1691) describe a linelike string that extends purely along the radial AdS 
coordinate, stretching upward from z = z m up to some turning point z = ((ir/2) and 
then returning back down to z = z m . The function ((a) will be smooth at the turning 
point only if <^ / (7r/2) = 0, which combined with (ITU]) and the requirement v x (tt/2) = 
implies that ((ir/2) = 1, i.e., the string turns around at the horizon. For this to occur, 
given an initial velocity profile v x (a), the value of A must be tuned in order for the 
numerical solution of (1701) to reach ( = 1 precisely at er = tt/2. So for these initial 
conditions, where Vq is a free parameter, A is not. Having determined C( CT )) one can 
proceed as before to the numerical integration of the equations of motion (1591) . A 
few representative quark trajectories are shown in Fig. [9j with v x (a) = t>ocos<7 and 
s — (1 — z) 1 / 2 (chosen to simplify the equation of motion and increase the stability 
of the numerical integration needed to obtain the initial profile ((&)). 



x 




0.2 0.4 0.6 0.8 1 1.2 1.4 



Figure 9: Quark trajectories for adjoint q-q configuration, for z m = 0.2 and v = v m 
(solid), vo = 0.75v m (dashed), v = 0.5v m (dashdotted). The corresponding energies 
are E/2M Test = 2.86, 1.38, 1.13. In this case no oscillating configurations are found. 

The existence of the two qualitatively distinct sets of initial conditions for the 
string describing the creation of a quark-antiquark pair has a direct field-theoretic 
interpretation: the product of a fundamental q and an antifundamental q can lead 
to a q-q pair either in the singlet or the adjoint representation of the SU(N C ) gauge 
group, and each of the above string configurations is dual to one of these. Indeed, the 
linelike string (169)) is precisely the system considered in jl5l [37J [62J, [18] to model a color 
source in the adjoint representation. Due to its extended nature, it sets up long-range 
supergravity fields that translate through the standard recipe of [2] into a long-range 
gluonic field profile, indicative of a source with net color charge. The completely 
pointlike string (1671) . on the other hand, sets up no long-range chromoelectromagnetic 
field and so describes the singlet. Given this correspondence, it is interesting that the 
AdS/CFT duality predicts that (at large N c and large A) the initial quark velocity is 
freely adjustable in the adjoint, but not the singlet, configuration. 
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4.3 Limiting velocity 

In the previous subsection we have learned that the initial velocity of the quark 
and antiquark at the moment of the pair's creation is bounded above by the speed 
v m defined in fl68|) . The reason for this is easy to understand on the string theory 
side of the duality, since, as we have noted above, the coordinate velocity v x = v rn 
corresponds to a proper velocity V x equal to that of light at the position of the string 
endpoints, z = z m [35J. The interesting feature is that it is v x , and not V x , that 
corresponds to the gauge theory velocity 

We would naturally expect the restriction to subluminal velocities, v < v m , to 
apply to more general string configurations. In complete analogy with the point 
particle case, the easiest way to deduce this restriction is to go back to the Nambu- 
Goto action (j2J), and observe that the requirement that it be real (i.e., that the 
string worldsheet be timelike) imposes a bound on physically realizable embeddings. 
Indeed, working for simplicity in the static gauge r = t, a — z, it is easy to see that 
the Nambu-Goto square root is real only as long as the embedding function X(z,t) 
satisfies 

< h -4nr — , (72) 



Of 



2 

M ff ) sin a 



where a = Z(dX/dt, dX/dz). For a string that moves and stretches along a single 
direction x, as we have considered up to now, a = or tt and this reduces to 

(t)-H(f) 2 )- 

It might seem peculiar that, as long as the string segment under consideration is 
not vertical (dX/dz ^ 0), the bound (1731) allows the proper velocity of the segment, 
(l/\/h)dX/dt, to exceed the speed of light (by an amount that becomes arbitrarily 
large in the limit \dX/dz\ — > oo). One should note, however, that this is a gauge- 
dependent statement, because unless the string segment is vertical, x is not entirely 
transverse to it, and the longitudinal component of the string motion is of course 
unphysical (in particular, motion along x is entirely unphysical in the limit \dX/ dz\ — > 
oo, where the string becomes horizontal). From 072)) we can see that, for purely 
transverse motion, a = ±7r/2 (as in the meson configurations we will consider in 
Section EJ), the proper velocity is indeed required to be less than unity, i.e., the 
coordinate (or, equivalently, gauge theory) velocity is bounded by \/I — z 4 . 

The argument of the preceding paragraph applies to a generic point in the interior 
of the string, but the situation is different for the endpoints, where motion along the 
body of the string is physical. As we know, a string endpoint on the D7-branes 
is dual to a quark (or antiquark), so the velocity v of the latter must necessarily 
respect the bound fl72|) . or, for collinear motion, f!73|) . Evaluating this last equation 
at z = z m and parametrizing as in Sections [2] and [3] the dXj ^-dependence in terms 
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of the momentum density II = 11^ given by ([3]), which controls the external force 
F = (a/A/27t)II applied to the quark, we deduce that v is bounded by 



h v 2 

v < 



y/(h m -ztIP)(l + z^P) ^(v 2 m - XF 2 /4n 2 m^)(l + XFyi^m^) ' 

(74) 

where h m = h(z m ) and [] denotes the expression within brackets in (1281) . 

For the case we have considered in the previous two subsections, where the string 
endpoint is free, corresponding to a quark that evolves only under the influence of 
the plasma, (174]) is precisely the statement that v < v m . As expected, we see that 
this bound (or, more generally, the analogous bound deduced from (1721 ) applies not 
only to quark-antiquark configurations, but also to isolated quarksjj 

For the case where the quark is externally forced, on the other hand, the bound 
( 1741) becomes less restrictive: as F increases, the quantity on the right grows mono- 
tonically, and in fact diverges when the force approaches the critical value _F cr i t = 
(\/ r \/2Tr)h m / mentioned already in Section [3721 (this is the same as the divergence 
seen in (172]) and (1731) when dX/dz — > oo). From this perspective alone, then, it would 
seem possible to take the quark to velocities larger than v m while exerting a force on 
it (even though, after release, one would again have v < v m ). 

It is not guaranteed, however, that this possibility is realized in practice. In the 
forced stationary case considered in [7J [8], for instance, it is found that v > v m would 
necessarily require F > F crit , and is consequently unattainable [20]. More generally, 
the question of whether v m is limiting or not is dynamical in nature, and essentially 
depends on the form of the thermal dispersion relation for the quark. In the discussion 
following (l5"2l we noted that, based on our zero-temperature results of Section [2J, the 
quark's intrinsic energy E q (and momentum p q ) should diverge as F — > F CTit . Given 
the connection between the 1 — v 2 factor in the denominator of (1521) and the Lorentz 
invariance of the metric (111) at T = 0, it is natural to expect it to be replaced by 
h m — v 2 = — v 2 at finite temperature, which would imply that E q — > oo (and 
p q — > oo) as v —* v m , meaning that v > v m is physically unattainable. The results we 
will obtain in Section 15.31 appear to support this expectation, but it would certainly 
be nice to be able to show this directlyjj 



4.4 Transition to asymptotic regime 

For the non-oscillatory string trajectories with the one-parameter family of initial 
conditions (1641) . it was found in [7J that the late-time quark motion is independent 
of the energy E of the configuration, and coincides with the behavior expected from 
a particle with relativistic dispersion relation p oc vj\J\ — v 2 , subject to a damping 

8 This observation has also been made very recently in |30j . Simultaneously, v m has been shown 
to emerge as a limiting velocity directly from the microscopic meson dispersion relation [39J. Both 
of these works appeared while the present paper was in preparation. 

9 For mesons this was done recently in |39j . which appeared while this paper was in preparation. 
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force 



dp 
dt 



VP , 



(75) 



with fi a ^-independent friction coefficient that was tabulated in [7J for various values 
of z m , and which is consistent with the drag force (C?T|) obtained in [7J, El Ej for the 
stationary configuration. 

The agreement between the analytic and late-time numeric results is most cleanly 
seen if instead of comparing graphs of x(t) or v (t) for the quark (where one would 
need to look at t — > oo), one examines the plots of t> (x)0 where the analytic behavior 
for constant /i takes the simple form 



which is linear with slope —\x near the final rest point x = Xqo (whose value is meant 
to be adjusted to fit the data). These plots are shown in Fig. [TUJ where it is seen that 
the late-time behavior is well-described by (T76]) also for oscillating trajectories and 
for the more general singlet configuration (l67j) . as well as for the adjoint configuration 
(169]) . Notice that the agreement holds in the late-time regime where (176]) reduces to a 
linear expression, but, as in the cases studied by the authors of [TJ, it extends beyond 
the range covered by their general quasi-normal mode analysis, because the velocity 
of the quark is not necessarily small. 

In Fig. [TUJ we see that there is an initial period where the behavior differs from 
the late-time frictional evolution (176]) . This difference is clearly more significant for 
the singlet than the adjoint case. The time that must elapse before the asymptotic 
behavior sets in becomes arbitrarily large for singlet configurations that are close to 
being oscillatory. In these cases, essentially all of the energy of the quark is lost 
not through the constant-/! frictional force due to the plasma, but as a result of the 
chromoelectromagnetic force exerted by the antiquark. 

On the string theory side of the duality, the issue is that, as the initially pointlike 
string grows and falls toward the black hole, it takes some time before it is close 
enough to z = 1 to be deformed into (two juxtaposed copies of) the asymptotic -i 
shape discussed in Section 13.31 This picture seems rather close to the phenomeno- 
logical discussion given in [50] (in the context of collisional energy loss): when the 
singlet quark- antiquark pair is formed within the plasma, there is a delay before the 
interaction between the newly created sources and the plasma can set up the long 
range gluonic field profile that is responsible for the late-time dissipation. 

To examine in more detail the transition to the late-time behavior (1751) . for a 
variety of trajectories we have determined the point (xf, Vf) beyond which the numeric 
v(x) curve agrees with the analytic curve (176|) to the accuracy indicated by the fraction 
/. Even though, judging by the effect of halving the grid spacing, our numerical 
results seem to be accurate to at least 1%, the precision with which we can determine 
(xf,Vf) is limited by the uncertainty in the fitting parameter Xqo, which we estimate 

10 We thank Antonio Garcia for suggesting this. 




(76) 



41 



V 




Figure 10: Quark evolution (velocity as a function of traveled distance, in units of 
1 / 7rT) for five different initial conditions. To explore the neighborhood of the charm 
quark, the mass parameter has been chosen as z m = 0.2, corresponding to a limiting 
velocity v m = 0.9992. The dotted curves show the results of our numerical integration, 
contrasted against fits in solid red that use the analytic expression (1761) . with the value 
\i = 0.25 obtained in [7J and an optimal choice of the stopping distance x^. The 
three dotted curves starting at the same point describe singlet configurations (and 
therefore have Vq = v m , as explained in the main text). The green, magenta and purple 
curves correspond respectively to total energies E/2M mst = 2.45, 2.45, 1.01 and initial 
string velocity profiles v cos a, v cos 3 a and v cos a (leading to Xoo = 4.15, 3.07, 0.12). 
Notice in particular that the purple curve describes a situation where the quark and 
antiquarks turn around and come to rest while approaching one another. The two 
remaining curves arise from adjoint configurations with different energies and initial 
quark velocities: the case E/2M rest = 1.12 and vq = 0.5v m (leading to Xoo = 2.07) is 
shown in dark blue, while E/2M rest = 1.05 and vq = 0.35v m (xqo = 1.4) is shown in 
light blue. 

to be of order 5-10%. It does not make much sense therefore to consider a value of / 
smaller than this. 

A representative sample of our results for / = 0.1 and / = 0.05, in the case 
of singlet configurations is shown in Fig. [TTJ The two sets of data have somewhat 
different functional forms, but are consistent with one another within the rather large 
margin of error. The general tendency is for xj to approach zero as v f — > v rn . 

For adjoint configurations, we find that the numeric v(x) curves are within 5- 
10% of (!76|) already at the start of the evolution, so, to be consistent with our pre- 
established criterion, in this case we must identify the transition length as essentially 
Xf — 0. It is still worth noting, however, that there is an initial period where the 
functional form of the numeric and analytic plots is different, as can be seen in 
Fig. [lOj This difference is negligible for small initial velocities, and becomes more 
pronounced (both in magnitude and in duration) as vq increases. The transition 
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length identifiable at the 3-4% level thus follows a trend opposite to the one for the 
singlet case. Nonetheless, up to the velocities Vf ~ 0.5 that we have been able to 



explore (corresponding to v 
Fig.QU 



it remains smaller than the singlet Xf shown in 
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Figure 11: Transition distance Xf (in units of 1/vrT) as a function of the transition 
velocity Vf, for / = 0.1 (red) and / = 0.05 (magenta). The size of the dots gives a 
rough indication of the margin of error. 



A question worth considering is whether the transition to the regime where the 
quark experiences a constant drag coefficient occurs right after the quark and an- 
tiquark are screened from each other by the plasma, or if there is an intermediate 
regime where the quark moves independently from the antiquark but nevertheless 
feels a drag force that differs from the stationary result of [TJ El [9] , as we found when 
applying an external force in Section 13.21 To answer this question, we will determine 
the relevant screening length in the next section. 

It would also be very interesting to develop an interpretation of energy loss for 
the evolving quark-antiquark pair in parallel with the picture for an isolated quark 
proposed by Mikhailov [52] and developed further in Sections I2ll3l An effort in this 
direction was in fact made in the work [53] , but regrettably we do not understand its 
use of trajectories that lie outside of the string worldsheet. We suspect that a treat- 
ment based on null curves on the worldsheet should be possible. Of course, progress 
is again hampered by the lack of an analytic solution describing the back-to-back q-q 
evolution considered in this section (or even its counterpart at zero temperature). 
Another potentially confusing issue is the fact that, in contrast with the isolated 
quark case, there are now two string endpoints at z — z m , and a priori it would 
be possible for null trajectories to 'bounce' repeatedly between them. One should 
however bear in mind that, at least beyond some finite interval of time, this would be 
prevented by the formation of a worldsheet horizon analogous to the one discussed in 
Sections CTO 



43 



5 Quark- Antiquark Potential 



5.1 Review of earlier results 

The potential E(L) for an infinitely massive static quark and antiquark in a strongly- 
coupled J\f = 4 SYM plasma was determined in using the dual description of 
the pair in terms a string with both of its endpoints on the AdS boundary z — 0. 
As expected, for small quark-antiquark separation L, the effects of the plasma are 
negligible and the potential matches the zero-temperature result [65] 

E(L) = -WSp . (77) 

As the separation grows, however, the effects of the plasma progressively screen the 
quark and antiquark from one another, and as a consequence raise the system's energy 
above the Coulombic behavior. The fl-shaped string embedding that is employed for 
the calculation exists only up to a maximal separation L max = 0.865/VT, and its en- 
ergy exceeds that of the disconnected solution beyond the somewhat smaller distance 

= 0.755/ttT. In the N c , A — > oo limit where the string-theoretic calculation is easy 
to carry out, the potential has a kink at L*, and beyond this point the string path 
integral is dominated instead by a configuration with high worldsheet curvature which 
for large L describes graviton exchange between two disconnected strings, leading to 
a tail E(L) T exp(-L/L gap ), with L gap = 0.428/vrT [38]. 

This computation can easily be extended to the case where the q-q pair moves 
with velocity v with respect to the plasma, to obtain the energy E(L,v) of the pair 
[28| [27]. For any v, the potential reduces to (1771) at small L. The main overall effect 
of increasing the velocity is to move the non-Coulombic portion of the E(L, v) curve 
down and to the left, leading in particular to the identification of the two screening 
lengths = L max for v > 0.447. 

The function L m ^(v) was determined numerically in [321 128] (and a related length 
was plotted in [31] )o Over the entire range < v < 1 its behavior may be approxi- 
mated as [28] 



W-)-^(l"- 2 ) 1/3 , (78) 
while in the ultra- relativistic limit, it can be shown analytically that [32] 

T / n 1 3^ 3 / 4 4tt 3 / 2 . 2U/4 0.743,., 2U/4 f 

^max(^) -> , , sn (1 - v 2 ) 1/4 ~ (1 - w 2 ) 1/4 for V -> 1. (79) 

V ; vrT T(l/4) 2 V 1 ttT v ; K ' 

The full curve L max (v) does not deviate far from this asymptotic form, so a decent 
approximation to it is obtained by replacing 0.743 — > 0.865 in (1501) . to reproduce the 



11 For the corresponding weak coupling calculation in QCD, see, e.g., |66j and references therein. 
12 The works [32l [27] additionally explored the dependence of L max on the angle 8 between the 
direction of motion and the dipole axis (which was fixed at 9 = tt/2 in [28 J, finding it to be weak. 
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correct value at v = (at the expense of introducing a 16% error as v — > 1) [32J: 

. , , 0.865,. 



,2U/4 



ttT 



50) 



A comparison between the two approximations (I7H1) and (|80|) is shown in Fig. | 
overall, the exponent 1/3 is better than 1/4 in the sense that it implies a smaller 
squared deviation from the numerical results, even though 1/4 leads to a smaller 
percentage error in the range v > 0.991 (7 > 7.3). An attempt to better parametrize 
the deviation away from the ultra-relativistic behavior was made in [36]. In any case, 
one should bear in mind that the region of principal interest at RHIC is not really 
v — > 1, but ryv ~ 1 (see, e.g., [6] 




Figure 12: Screening length L max (in units of l/27rT) as a function of velocity (in 
black) compared against the approximations fITBl (in red) and flHUj) (in blue). 



It is interesting to consider how these results generalize to the case where the quark 
and antiquark have a large but finite Lagrangian mass m, which as reviewed in Section 
[2] corresponds to letting the flavor D7-branes extend up to a position z m > given 
by (l28|) . As explained in Section [3T2l from the phenomenological perspective we are 
mostly interested in values of the mass parameter z m in the range z m /zh ~ 0.2 — 0.4, 
which corresponds to the charm mass. 



5.2 Finite mass at zero temperature 

In a static configuration, the equation of motion for the string amounts to the state- 
ment that the momentum density IT^. given by ([3]) is independent of x. In the pure 
AdS geometry (i.e., in the metric ([T]) with T = 0) the resulting embedding is 



/"2max/2 J/ 1 

X(z) = ±z max / p , (81) 

where 2 max e [z m , oo) denotes the turning point of the fl-shaped string and the plus 
or minus sign applies respectively to its right and left half. From this it follows in 
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particular that the distance between the string endpoints at z = z n 
separation) is 



2z„ 



^max / Zm. 



c 2 Vc 4 



2z 
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; 4'2 



(i.e., the q-q 



(82) 

where B denotes the incomplete Euler beta function. The energy of the fl-shaped 
string, renormalized by subtracting the energy of two purely radial strings, follows 
from ([2]) as 
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(83) 



The corresponding results for the case of a moving q-q pair are of course determined 
in terms of these by Lorentz invariance. 

Using (182]) and (1531 . it is easy to see that for L >> \f\jm (i.e., z max S> z m ) the 
quark-antiquark potential E(L) approaches the same Coulombic form (177)1 as in the 
infinitely massive case (z m = 0). On the other hand, in the short distance limit L — > 
we have z max — > z m and so the potential is linear, 



E(L) = 2m 



it m I ( mV 
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These results are consistent with the form of the gluonic field around an isolated 
finite-mass quark 
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(where Oi?2 = Tr{F 2 + . . .} /Agy M denotes the operator dual to the dilaton field [68]). 
which is Coulombic at long distances and non-singular at the location of the quark. 

In both (1841) and (|85l) we are seeing the effect of a color charge distribution that 
cloaks the fundamental sources, or in other words, of the gluon (and scalar, etc.) cloud 
that surrounds the would-be bare quark to turn it into a 'dressed' or 'composite' 
quark, in line with the inferences we made in Section I2.2L Notice from ( 1851) that 
the size of this cloud is not the Compton wavelength of the quark, 1/m (indicating 
that it is not made of virtual quarks and antiquarks), but the much larger scale 
\/A/m, which is the characteristic size of the deeply bound microscopic mesons of the 
theory [691 EH E9] — so perhaps it would be more appropriate to refer to this charge 
distribution as a 'meson cloud.' 
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5.3 Finite mass at finite temperature 



We will consider first the case where the quark- ant iquark pair is static with respect 
to the plasma. Repeating the calculation of the previous subsection for T > 0, the 
separation between the quark and antiquark is found to be 



L = 2z„ 



2max / Zm 



d( 



<r 4 



with ( h 
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J Zh, and the renormalized energy of the pair is given by 



(86) 
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From (1861) and (1871) . the short-distance behavior of the potential is now 
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with z m related to the quark Lagrangian mass m through (128]) . By comparing with 
numerical plots (as in Fig. [13] below), this linear expression can be seen to give a good 
approximation of the actual potential at least up to L ~ z m /2. 

The constant term in (I88p is clearly nothing but (minus) the thermal rest mass 
([55]) of the isolated quark and antiquark, which was subtracted in our choice of renor- 
malization. For any L, the force F pot (L) = —dE/dL deduced from the potential 
balances the external force 



dp _V\ 
t[L} ~ dt ~ 27 




59) 



required to hold the string endpoint in place, and the coefficient of the linear term in 
f ]88]) correctly encodes the L — > (z max — ► z m ) limit of this force. Using (JT]) and (168]) . 
this can in turn be rewritten in the form 



7l 

2 



(90) 



which coincides with the drag force (I37p computed with a trailing string in [7J [SJ, 
evaluated at the velocity v = v m . This agreement has a simple interpretation. As 
explained in [20], the drag force calculation of [7] [S] has a limited range of validity 
( z m — z h(^ — v 2 ) l / A for a given v, or equivalently, v <v m for a given z m ), arising from 
the condition that the electric field F 0x needed to keep the string endpoint moving at 
constant speed does not exceed its critical value at z m , Fq"* = (a/A / '2tt) h(z m ) (zh / z m ) 2 . 
This critical electric field is thus equivalent to the maximal attainable drag force, 
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F drag (f m ). The fact that it agrees with the right-hand side of (!90|) . then, is not 
surprising, because, by definition, F ext (0) encodes the force needed to pull apart the 
endpoints of a zero-size string at z = z m , i.e., the value of F^} 1 . In gauge theory- 
language, this is the force required to nucleate a quark- ant iquark pair. 

For L > z m /2, the potential is screened and deviates significantly from the linear 
form flHHj) . Its behavior is shown in Fig. [131 As can be seen there, the location 
of the screening lengths L* and L max shift to the left with increasing z m (e.g., for 
z m /z h = 0,0.2,0.4,0.75 we respectively find vrTL, = 0.755,0.746,0.707,0.511 and 
7rTL max = 0.865, 0.861, 0.828, 0.616). 
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Figure 13: Quark-antiquark potential (in units of T^/X/A) in the static finite-mass 
case, for z m /z h = 0.001 (red), 0.2 (light green), 0.3 (orange), 0.4 (light blue), 0.5 
(blue) and 0.75 (purple). The separation L is given in units of l/2irT. See the main 
text for discussion. 



Let us now go on to consider the case where the quark-antiquark pair moves 
through the plasma. For simplicity, as in [28] we will restrict ourselves to the case 
where the motion is perpendicular to the dipole axis (other angles have been studied 
in [321 ESI EZ])- Orienting the former along x and the latter along y, the expression 
for the q-q separation can be copied directly from Eq. (37) of [25] . It is given by 

f r hm dh 

L ^ v) = Sr 1 I (91) 

27rT A min (i _ h)L (h - v 2 )h - (l - h)hP 



where 

„2ttz 



f y = <n 



y ft, /VT2 



is a rescaled version of the force F y needed to hold the quark and antiquark in place, 
h(z) is the function appearing in the metric ([1]), 



and we have taken into account the finiteness of the quark Lagrangian mass m by 
changing the upper limit of the integration range from 1 to h m = h(z m ). The form of 
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expression (f9~TT) reflects the fact that, in spite of its motion with respect to the black 
hole, the string remains completely vertical, in contrast with the trailing configuration 
studied in [H [SJ. In the gauge theory, this translates into the interesting property 
that, unlike isolated quarks, mesons feel no drag [3T | |32 }J 28]. 

Similarly, the energy of the q-q pair in its rest framd 13 ! can be read from Eq. (38) 
of [28] as 



E(fy,v) 



dh(h — 1> 2 )7 



(1 - h)^(h-v 2 )h- (] 



h)hfl 



dh 



[i-hy. 



(93) 



The second term in this expression represents the energy of two disjoint strings trailing 
their boundary endpoints as in [7J [8] , dual to an unbound quark and antiquark (which 
do experience a drag force). This contribution is subtracted in order for E = to 
correspond to the energy beyond which the system would become unbound. 
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-2.5 L ' ' ' -2.5 L *» ' 

Figure 14: Quark-antiquark energy (in units of Ty/X/A) function of separation 
(in units of 1/2ttT), for (a) fixed v = 0.45 and z m = (black), 0.2 (red), 0.25 (orange), 
0.5 (light green), and 0.75 (green); (b) fixed z m = 0.2 and v = (black), 0.2 (red), 0.4 
(orange), 0.6 (light green), 0.8 (green) and 0.9 (bright green). See text for discussion. 

Carrying out the integrals in (1911) and (1931) numerically, one obtains plots for the 
potential such as those shown in Fig. [TH For small separation the behavior is again 
linear, 



E(L) = 




The two terms in this expression have the same interpretation as was given for (1881) 
in the paragraph surrounding (1591 - (|9~01) . 

From Fig. [TH one also sees that the value of z m has an effect on the location of 
the screening lengths L*, L max (and the velocity beyond which L* = L max ). The 

13 For a discussion of subtleties in the calculation of the pair energy in the rest frame of the plasma, 
see EH1- 
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behavior at fixed v and varying z m shown in Fig. [T4"r is qualitatively the same as 
we saw already for the static configuration in Fig. [TSl The dependence at fixed z m 
and varying v, which is more relevant for the phenomenological situation, is shown 
in Fig. [lib. 




v 



Figure 15: (a) Screening length L max (in units of 2ttT) as a function of velocity for 
z m / %h = 0.2 (in black) compared against the z m > fits (|95|) (in red) and (196]) (in light 
blue), (b) Expanded version of the same plot, showing that (|95|) gives a relatively 
good approximation up to velocities very close to v m , where the asymptotic linear 
behavior (1971) (dotted dark blue) sets in. 

In Figs. [T5l and [T6| we see that when z m /zh = 0.2, and to a lesser extent, when 
Zm/zh = 0.4, the screening length L max (v) is still relatively well approximated in the 
full range < v < v m by the natural modification of the z m = fit (|78|) . 

W-)-^(^-- 2 ) 1/3 - (95) 

This approximation becomes worse as z m /zh is further increased. (In all cases, the 
fit analogous to flHOl . 

W^-^K.-- 2 ) 1/4 , (96) 

does a poorer job than (j95l) .) 

An important difference with respect to the case where the quarks are infinitely 
massive is that now the 'ultra-relativistic' region would refer to the limit where the 
pair velocity approaches the limiting velocity v m < 1 given in fl68|) and discussed 
further in Section 14.31 From this alone it is clear that the asymptotic formula flHOl) 
no longer holds. The behavior in this region can still be determined analytically, and 
turns out to be 

r J_ 1 Zj[l - (Z m/ 'Z h f - J\ 

L — 27^(1-^)3/4 2nT _ {Zm , Zh y ' ^ 3<) 

where we see that the 1/4 exponent in ( |79l) changes to 1 at finite quark mass. As 
shown in Figs. fT5b and fT6b. for values of z m /zh in the neighborhood of the charm 
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Figure 16: (a) Screening length L max (in units of 2tiT) as a function of velocity for 
Zm/zh = 0.4 (in black) compared against the z m = fits (178"|) (in red) and (1H0|) (in 
light blue), (b) Expanded version of the same plot, showing that the region in the 
neighborhood of v m where the asymptotic linear behavior (l97j) (dotted dark blue) is 
a significantly better approximation than (1781) has grown larger in comparison to the 
Zm/zh = 0.2 case depicted in Fig. [151 

mass, the functional form ( 1971) applies only for velocities extremely close to the v m 
endpoint. As z m /zh is further increased, however, there is a tendency for ( 1971) to 
apply to a larger fraction of the < v < v m interval. 

In [32| [27] , the result (IHUl) for the screening length in the infinitely massive case 
was turned into a tentative prediction for the dissociation temperature of charmonium 
or bottomonium traversing the quark-gluon plasma, Tdi SS oc (1 — v 2 ) 1 ^. (As reviewed 
in Section 147X1 for < v < 0.991 the exponent implied by the AdS/CFT calculation 
is really closer to 1/3 than to 1/4.) The argument for this velocity scaling is heuristic, 
and relies on comparing the screening length against the natural size of the bound 
state in question (see, e.g., [71] and references therein). Applying the same logic to 
the fits (1951) or (l9"7j) for the screening length in the case of masses in the neighborhood 
of the charm quark {z m jz h ~ 0.2-0.4, as in Figs. [T5TfT6l) . one would infer that T diss oc 
(v^ — v 2 ) n , with n ~ 1/3 for general v, and n — 1 in the v — > v m limit (the latter 
being the preferred fit only in the restricted range v > 0.998 or v > 0.98, depending 
on the value of z m ). 

While this paper was in preparation, the work [39] appeared, which includes an 
explicit derivation of the dispersion relation for mesons, described as open string 
modes on the D7-branes, following [69]. For large meson momentum p, their calcula- 
tion resulted in E = v m p, with v m the limiting velocity defined in (!68l) and discussed 
further in Section 14.31 of the present paper (as well as in [30] , which also appeared 
while our work was being written up). The authors of [39] emphasized that the z^- 
and 2 m -dependence of v m implies a bound on the temperature at which a meson with 
velocity v can exist within the plasma. Interpreting the maximal allowed temperature 
as a dissociation temperature, they found the scaling Tdi SS oc (1 — v 2 ) 1 ^, just like the 
heuristic argument based on the potential for infinitely-massive and ultra-relativistic 
quarks had suggested [321 [27J . 

Given the results of [2B] and the present subsection, we are not quite sure what 
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to make of this agreement. As explained above, for finite mass the quark- ant iquark 
potential does not scale as (1 — v 2 ) 1 I A even in the 'ultra-relativistic' region (which is 
now v —> v m ). So the fact that the heuristic argument in [321 [21] based on the potential 
suggests the same scaling with velocity as the direct calculation in [27] might be a 
coincidence, or it might indicate that there is a direct physical connection between 
the dissociation temperature and the potential only when the latter is computed for 
infinitely- massive quarks. In any case, it should be emphasized that the 1/4 scaling 
follows not from the detailed computation of the dispersion relation in [27], but from 
the temperature-dependence of the limiting velocity (|68|) . 

5.4 Screening length vs. transition distance 

In Fig. [TT] at the end of the previous Section, we determined the location (xf,Vf) 
beyond which a quark and antiquark created (with zero total momentum) within 
the plasma begin to slow down at the asymptotic rate specified by the constant 
friction coefficient obtained for an isolated quark in [7J [8]. For q-q separations larger 
than 2xj, then, the two particles are definitely oblivious to one another and evolve 
independently. It is therefore interesting to compare the transition distance Xf against 
(half of) the screening length L max determined in the present section, which gives a 
measure of the minimal distance at which the quark and antiquark can decouple. 

For singlet configurations, the result of this comparison is shown in Fig. [TTJ where 
we see that the two separations are of comparable magnitude and scale with velocity in 
a similar manner0 Notice that this is in spite of the fact that the two relevant string 
configurations are quite different, with motion respectively along and perpendicular 
to the plane in which the string extends. We conclude then that the transition 
to the constant-drag-coefficient regime takes place immediately after the quark and 
antiquark lose contact with one another. That is to say, unlike what we found for 
the forced isolated quark in Section 13.21 here there is no intermediate stage where the 
quark and antiquark decelerate independently from one another at a rate that differs 
substantially from the asymptotic result of [7J1H]. 

The initial stage, where as seen in Fig.[10]the quark and antiquark dissipate energy 
at a rate larger than (1371) . corresponds to the period when the string has not fallen 
close enough to the black hole horizon. In other words, to a first approximation, the 
quark and antiquark evolve in this region as they would in the absence of the thermal 
plasma. To the extent that this result, with all its simplifying assumptions, might 
conceivably be extrapolated to the experimental context, this initial stage would not 
differ significantly between heavy ion and proton-proton collisions. 

For adjoint configurations, the quark- ant iquark potential is expected to be sup- 
pressed by a factor of 1/N 2 relative to the singlet case (see, e.g., [72] and references 
therein). On the string theory side of the duality, this is visible in the fact that the 
'line' initial conditions (1691) can be regarded as describing two disjoint strings extend- 

14 A similar agreement was reported very recently in [40] . which appeared while this paper was in 
preparation. 
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Figure 17: (Half of) Screening length L max (v) (green) compared against the transition 
distance Xf(v) defined at the end of Section 0], with z m /zh = 0.2, for / = 0.1 (red) 
and / = 0.05 (magenta). The vertical axis is in units of 1/nT. 



ing all the way from the D7-branes to the horizon. In the classical limit where our 
calculations are carried out, these two halves do not influence one another, because 
signals propagating between them would have to pass through the horizon, and their 
splitting/joining/annihilating is suppressed by powers of the string coupling g s , which 
is taken to zero. At this level of approximation, then, the quark and antiquark do 
not see one another, so the relevant 'screening length' is L s = 0. This is consistent 
with the result Xf = we obtained for the adjoint case in Section WM 

The fact that in the adjoint case the quark and antiquark evolve independently 
from the start makes this case directly comparable to our study of isolated quarks 
in Section 13.21 Interestingly, the corresponding results seem to point in opposite 
directions: whereas for an initially static quark that is externally forced we identified 
a short period after release where the quark loses energy at a rate significantly smaller 
than the stationary/asymptotic rate (1371) . for a quark that is created with a certain 
initial velocity and gluonic field profile and thereafter moves only under the influence 
of the plasma we found an initial rate of dissipation that is somewhat larger than 



It seems, then, that the initial conditions can, to a certain extent, influence the 
subsequent evolution, which does not seem altogether surprising given the nonlinear 
nature of the medium. The application of an external force in Sect ion 13721 clearly has a 
significant effect on the shape of the string at the time of release, or in other words, on 
the gluonic field configuration surrounding the quark. In view of our results for quark- 
antiquark evolution in this and the previous section, the initial period over which this 
disturbed quark behaves unlike what one would expect based on [7J El [9] should most 
likely be interpreted as reflecting the time it takes the quark to stabilize the gluonic 
fields around it, or, equivalently, to move far enough away from the disturbed region 
to be effectively screened from it. In any case, one should of course not lose sight of 
the fact that the actual experimental situation resembles the setup of Section H] more 
closely than that of Section 13.21 



fTJ). 



53 



Acknowledgements 



It is a pleasure to thank Alejandro Ayala, Elena Caceres, Jorge Casalderrey-Solana, 
Antonio Garcia and Guy Paic for valuable discussions. We are also grateful to Antonio 
Garcia for collaboration in the initial stages of this project. We additionally thank 
Karen Ramirez for help with Figs. Q] and [71 This work was partially supported by 
Mexico's National Council of Science and Technology (CONACyT) grants 50-1551 and 
SEP-2004-C01-47211, as well as by DGAPA-UNAM grants IN116408 and IN104503. 



References 

[1] J. M. Maldacena, "The large N limit of superconformal field theories and super- 
gravity,' 1 ^A^y_^Theor^fath_Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 
(1999)] |arXiv:hep-th/9711200| . 

[2] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, "Gauge theory cor- 
relators from non-critical string theory," Phys. Lett. B 428, 105 (1998) 
|arXiv:hep-th/980210"9] ; 

E. Witten, "Anti-de Sitter space and holography," Adv. Theor. Math. Phys. 2, 
253 (1998) |arXiv:hep-th/9802150l . 

[3] O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri and Y. Oz, "Large 
N field theories, string theory and gravity," Phys. Rept. 323, 183 (2000) 
|arXiv:hep-th/9905 lTT] . 

[4] W. A. Zajc, "The Fluid Nature of Quark-Gluon Plasma," larXiv:0802.3552l 
[nucl-ex] ; 

B. Muller, "From Quark-Gluon Plasma to the Perfect Liquid," Acta Phys. 
Polon. B 38, 3705 (2007) |arXiv:0710.3366l [nucl-th]]; 

E. Shuryak, "Emerging theory of strongly coupled quark-gluon plasma," 
arXiv:hep-ph/0703208| 

D. d'Enterria, "Qua rk-gluon matter," J. Phys. G 34 (2007) S53 
|arXiv:nucl-ex/061l6~12] ; 

C. A. Salgado, "Heavy ions theory review," Acta Phys. Polon. B 38 (2007) 975 
arXiv:hep-ph/0609i72] ; 

E. V. Shuryak, "Strongly coupled quark-gluon plasma: The status report," 
arXiv:hep-ph/0608177| 

M. J. Tannenbaum, "Recent results in relativistic heavy ion collisions: From 'a 
new state of matter' to 'the perfect fluid'," Rept. Prog. Phys. 69 (2006) 2005 
arXiv:nucl-ex/060300~3] ; 

B. Muller and J. L. Nagle, "Results from the Relativistic Heavy Ion Collider," 
|arXiv:nucl-th/06020~29| 

M. Gyulassy and L. McLerran, "New forms of QCD matter discovered at 
RHIC," Nucl. Phys. A 750, 30 (2005) |arXiv:nucl-th/0405013] . 



54 



[5] K. Adcox et al. [PHENIX Collaboration], "Formation of dense partonic matter 
in relativistic nucleus nucleus collisions at RHIC: Experimental evaluation by the 
PHENIX collaboration," Nucl. Phys. A 757, 184 (2005) |arXiv:nucl-ex/0410003| ; 
I. Arsene et al. [BRAHMS Collaboration], "Quark gluon plasma and color glass 
condensate at RHIC? The perspective from the BRAHMS experiment," Nucl. 
Phys. A 757 (2005) 1 |arXiv:nucl-ex/0410020] ; 

B. B. Back et al, "The PHOBOS perspective on discoveries at RHIC," Nucl. 
Phys. A 757 (2005) 28 |arXiv:nucl-ex/0410022| ; 

J. Adams et al. [STAR Collaboration], "Experimental and theoretical challenges 
in the search for the quark gluon plasma: The STAR collaboration's critical 
assessment of the evidence from RHIC collisions," Nucl. Phys. A 757, 102 (2005) 
|arXiv:nucl-ex/0501009] . 

[6] F. Carminati et al. [ALICE Collaboration], "ALICE: Physics performance re- 
port, volume I," J. Phys. G 30 (2004) 1517; 

A. Morsch [ALICE Collaboration], "Jet quenching studies with the ALICE de- 
tector," Czech. J. Phys. 55 (2005) B333; 

L. I. Sarycheva, "Heavy ion collisions at CMS and quark gluon plasma," Czech. 
J. Phys. 52 (2002) C171; 

D. d'Enterria et al., "CMS physics technical design report: Addendum on high 
density QCD with heavy ions," J. Phys. G 34 (2007) 2307. 

[7] C. P. Herzog, A. Karch, P. Kovtun, C. Kozcaz and L. G. Yaffe, "Energy loss of a 
heavy quark moving through M = 4 supersymmetric Yang-Mills plasma," JHEP 
0607 (2006) 013 [arXiv:hep-th/0605158] . 



[8] S. S. Gubser, "Drag force in AdS/CFT," Phys. Rev. D 74 (2006) 126005 
[arXiv:hep-th/0605182] . 

[9] J. Casalderrey-Solana and D. Teaney, "Heavy quark diffusion in strongly coupled 



M = 4 Yang Mills," Phys. Rev. D 74 (2006) 085012 |arXiv:hep-ph/0605199 



[10] H. Liu, K. Rajagopal and U. A. Wiedemann, "Calculating the jet 
quenching parameter from AdS/CFT," Phys. Rev. Lett. 97 (2006) 182301 
|arXiv:hep-ph/0605178] . 

[11] G. Policastro, D. T. Son and A. O. Starinets, "The shear viscosity of strongly 
coupled N = 4 supersymmetric Yang-Mills plasma," Phys. Rev. Lett. 87 (2001) 
081601 |arXiv:hep-th/0104066] ; 



A. Buchel and J. T. Liu, "Universality of the shear viscosity in supergravity," 



Phys. Rev. Lett. 93, 090602 (2004) [arXiv:hep-th/0311175| ; 
P. Kovtun, D. T. Son and A. O. Starinets, "Viscosity in strongly interacting 
quantum field theories from black hole physics," Phys. Rev. Lett. 94 (2005) 
111601 [arXiv:hep-th/0405231] ; 

A. Buchel, J. T. Liu and A. O. Starinets, "Coupling constant dependence of the 



55 



shear viscosity in M = 4 supersymmetric Yang-Mills theory," Nucl. Phys. B 707, 
56 (2005) [arXiv:hep-th/0406264]. 



[12] A. O. Starinets, "Transport coefficients of strongly coupled gauge theories: In- 
sights from string theory," Eur. Phys. J. A 29 (2006) 77 [ arXiv:nucl-th/0511073| ; 
D. T. Son and A. O. Starinets, "Viscosity, Black Holes, and Quantum Field 
Theory," Ann. Rev. Nucl. Part. Sci. 57 (2007) 95 [arXiv: 0704 .02401 [hep-th]]; 
R. B. Peschanski, "Quark-Gluon Plasma/Black Hole duality from 
Gauge/Gravity Correspondence," larXiv:0710.0756l [hep-ph]. 

[13] A. Majumder, B. Muller and X. N. Wang, "Small Shear Viscosity of a Quark- 
Gluon Plasma Implies S trong Jet Quenching," Phys. Rev. Lett. 99 (2007) 192301 
|arXiv:hep-ph/0703d82| . 

[14] R. Baier, D. Schiff and B. G. Zakharov, "Energy loss in perturbative QCD," 
Ann. Rev. Nucl. Part. Sci. 50 (2000) 37 [arXiv:hep-ph/0002198] ; 
R. Baier, "Jet quenching," Nucl. Phys. A 715, 209 (2003) 



arXiv:hep-ph/0209038|; 



A. Kovner and U. A. Wiedemann, "Gluon radiation and parton energy loss," 
|arXiv:hep-ph/030415lj 

[15] J. Casalderrey-Solana and C. A. Salgado, "Introductory lectures on jet quenching 
in heavy ion collisions," larXiv: 0712 .34431 [hep-ph]; 

D. E. Kharzeev, "Theoretical issues in J/psi suppression," J. Phys. G 34 (2007) 
S445 [arXiv:hep-ph/0703259] ; 

R. Rapp and H. van Hees, "Heavy Quark Diffusion as a Probe of the Quark- 
Gluon Plasma," larXiv:0803.0901l [hep-ph]. 

[16] C. P. Herzog, "Energy loss of heavy quarks from asymptotically AdS geome- 
tries," JHEP 0609 (2006) 032 |arXiv:hep-th/0605191| ; 

E. Caceres and A. Giiijosa, "Drag force in a charged M = 4 SYM plasma," 
JHEP 0611 (2006) 077 [arXiv:hep-th/0605235] ; 

S. J. Sin and I. Zahed, "Ampere's law and energy loss in AdS/CFT duality," 
Phys. Lett. B 648 (2007) 318 [arXiv:hep-ph/0606049| ; 

T. Matsuo, D. Tomino and W. Y. Wen, "Drag force in SYM plasma with B 
field from AdS/CFT," JHEP 0610 (2006) 055 |arXiv:hep-th/0607178| ; 
P. Talavera, "Drag force in a string model dual to large-iV QCD," JHEP 0701 
(2007) 086 |arXiv:hep-th/0610179| ; 

E. Antonyan, "Friction coefficient for quarks in supergravity duals," 



arXiv:hep-th/0611235 



D. Berenstein and H. J. Chung, "Aspects of open strings in Rindler Space," 
larXiv:0705.3lT0l [hep-th]; 

A. Karch and A. O'Bannon, "Metallic AdS/CFT," JHEP 0709 (2007) 024 
|arXiv:0705.3870l [hep-th]]; 

C. P. Herzog and A. Vuorinen, "Spinning Dragging Strings," JHEP 0710 (2007) 



56 



087 |arXiv:0708.0609l [hep-th]]; 

A. O'Bannon, "Hall Conductivity of Flavor Fields from AdS/CFT," Phys. Rev. 
D 76 (2007) 086007 |arXiv:0708. 19941 [hep-th]]; 

K. B. Fadafan, ll R 2 curvature-squared corrections on drag force," 
larXiv:0803.2777l [hep-th]; 

J. F. Vazquez-Poritz, "Drag force at finite 't Hooft coupling from AdS/CFT," 
to appear. 

[17] E. Caceres and A. Giiijosa, "On drag forces and jet quenching in strongly cou- 
pled plasmas," JHEP 0612 (2006) 068 |arXiv:hep-th/0606134| ; 
E. Nakano, S. Teraguchi and W. Y. Wen, "Drag Force, Jet Quenching, and 
AdS/QCD," Phys. Rev. D 75 (2007) 085016 |arXiv:hep-ph/0608274| ; 
G. Bertoldi, F. Bigazzi, A. L. Cotrone and J. D. Edelstein, "Holography 
and Unquenched Qua rk-Gluon Plasmas," Phys. Rev. D 76 (2007) 065007 
arXiv:hep-th/0702225] . 



[18 
[19 
[20 

[21 

[22; 
[23 
[24 



M. Chernicoff and A. Giiijosa, "Energy loss of gluons, baryons and fc-quarks in 
an M = 4 SYM plasma," JHEP 0702 (2007) 084 |arXiv:hep-th/0611155| . 

S. S. Gubser, "Momentum fluctuations of heavy quarks in the gauge-string du- 
ality," Nucl. Phys. B 790 (2008) 175 [arXiv:hep-th/0612143] . 

J. Casalderrey-Solana and D. Teaney, "Transverse momentum broadening of 
a fast quark in a N = 4 Yang Mills plasma," JHEP 0704 (2007) 039 
|arXiv:hep-th/0701123| . 

V. Mazu and J. Sonnenschein, "Non critical holographic models of the thermal 
phases of QCD," larXiv: 071 1.42731 [hep-th]. 

C. Athanasiou, H. Liu and K. Rajagopal, "Velocity Dependence of Baryon 
Screening in a Hot Strongly Coupled Plasma," larXiv:0801.1117l [hep-th]. 

S. S. Gubser, D. R. Gulotta, S. S. Pufu and F. D. Rocha, "Gluon energy loss in 
the gauge-string duality," larXiv:0803.1470l [hep-th]. 

P. M. Chesler and A. Vuorinen, "Heavy flavor diffusion in weakly coupled M = 4 
super Yang-Mills theory," JHEP 0611 (2006) 037 [arXiv:hep-ph/0607148| ; 
S. C. Huot, S. Jeon and G. D. Moore, "Shear viscosity in weakly coupled N = 4 
super Yang-Mills theor y compared to QCD," Phys. Rev. Lett. 98 (2007) 172303 
|arXiv:hep-ph/0608062| ; 

J. P. Blaizot, E. Iancu, U. Kraemmer and A. Rebhan, "Hard-thermal-loop 
entropy of supersymmetric Yang-Mills theories," JHEP 0706 (2007) 035 



arXiv:hep-ph/0611393|; 



W. A. Horowitz and M. Gyulassy, "Testing AdS/CFT Deviations from pQCD 
Heavy Quark Energy Loss with Pb+Pb at LHC," l larXiv: 0706 .23361 [nucl-th]; 
W. A. Horowitz, "Ratio of charm to bottom Raa as a test of pQCD vs. AdS/CFT 



57 



energy loss," larXiv:0710.0595l [hep-ph]; 

W. A. Horowitz, "pQCD vs. AdS/CFT Tested by Heavy Quark Energy Loss," 
larXiv:0710.0703l [nucl-th]; 

S. Caron-Huot and G. D. Moore, "Heavy quark diffusion in QCD and M = 4 
SYM at next-to-leading order," JHEP 0802 (2008) 081 |arXiv:0801.2T73l [hep- 
ph]]. 

J. J. Friess, S. S. Gubser and G. Michalogiorgakis, "Dissipation from a heavy 
quark moving through J\f = 4 super- Yang-Mills plasma," JHEP 0609 (2006) 072 
|arXiv:hep-th/0605292| ; 

Y. h. Gao, W. s. Xu and D. f. Zeng, "Wake of color fields in charged M = 4 



SYM plasmas," arXiv:hep-th/0606266; 



J. J. Friess, S. S. Gubser, G. Michalogiorgakis and S. S. Pufu, "The stress tensor 
of a quark moving through M = 4 thermal plasma," Phys. Rev. D 75 (2007) 



106003 |arXiv:hep-th/0607022|; 



A. Yarom, "The high momentum behavior of a quark wake," Phys. Rev. D 75 



(2007) 125010 arXiv:hep-th/0702164|; 



S. S. Gubser and S. S. Pufu, "Master field treatment of metric pertur- 
bations sourced by the trailing string," Nucl. Phys. B 790 (2008) 42 
|arXiv:hep-th/0703090| ; 

A. Yarom, "On the energy deposited by a quark moving in an N=4 SYM plasma," 
Phys. Rev. D 75 (2007) 105023 |arXiv:hep-th/0703095| ; 

S. S. Gubser, S. S. Pufu and A. Yarom, "Energy disturbances due to a moving 
quark from gauge-string duality," JHEP 0709 (2007) 108 |arXiv:0706.0"2l3l [hep- 
th]]; 

P. M. Chesler and L. G. Yaffe, "The wake of a quark moving through a strongly- 
coupled M = 4 supersymmetric Yang-Mills plasma," Phys. Rev. Lett. 99 (2007) 
152001 [arXiv:0706.0368l [hep-th]]; 

S. S. Gubser, S. S. Pufu and A. Yarom, "Sonic booms and diffusion wakes gen- 
erated by a heavy quark in thermal AdS/CFT," Phys. Rev. Lett. 100 (2008) 
012301 [ arXiv:0706.4307l [hep-th]]: 

S. S. Gubser and A. Yarom, "Universality of the diffusion wake in the gauge- 
string duality," larXiv:0709. 10891 [hep-th]; 

S. S. Gubser, S. S. Pufu and A. Yarom, "Shock waves from heavy-quark mesons 
in AdS/CFT," larXiv:0711.1415l [hep-th]; 

P. M. Chesler and L. G. Yaffe, "The stress-energy tensor of a quark moving 
through a strongly-coupled N=4 supersymmetric Yang-Mills plasma: comparing 
hydrodynamics and AdS/CFT," larXiv:0712.0050l [hep-th]; 

J. Noronha, G. Torrieri and M. Gyulassy, "Near Zone Navier-Stokes Analysis 
of Heavy Quark Jet Quenching in an M = 4 SYM Plasma," larXiv: 0712 .10531 
[hep-ph]; 

S. S. Gubser and A. Yarom, "Linearized hydrodynamics from probe-sources in 
the gauge-string duality," larXiv:0803.0081 l [hep-th]. 



58 



[26] A. Buchel, "On jet quenching parameters in strongly coupled non-conformal 
gauge theories," Phys. Rev. D 74 (2006) 046006 [arXiv:hep-th/0605178| ; 
J. F. Vazquez-Poritz, "Enhancing the jet quenching parameter from marginal 
deformations," jarXiv:hep-th/0605296; 

F. L. Lin and T. Matsuo, "Jet quenching parameter in medium with chemical 
potential from AdS/CFT," Phys. Lett. B 641 (2006) 45 |arXiv:hep-th/0606136| ; 
S. D. Avramis and K. Sfetsos, "Supergravity and the jet quenching pa- 
rameter in the presence of R-charge densities," JHEP 0701 (2007) 065 
|arXiv:hep-th/0606190| ; 

N. Armesto, J. D. Edelstein and J. Mas, "Jet quenching at finite 't Hooft 
coupling and chemical potential from AdS/CFT," JHEP 0609 (2006) 039 
|arXiv:hep-ph/0606245l ; 

Y. h. Gao, W. s. Xu and D. f. Zeng, "Jet quenching parameters of Sakai-Sugimoto 
model," |arXiv:hep-th/0611217[ 

A. L. Cotrone, J. M. Pons and P. Talavera, "Notes on a SQCD-like plasma dual 
and holographic renormalization," JHEP 0711 (2007) 034 jarXiv:0706.2766l [hep- 
th]]; 

Y. h. Gao, W. s. Xu and D. f. Zeng, "Viscosity and jet quenching from holo- 
graphic model," larXiv:0707.0817l [hep-ph]. 

[27] H. Liu, K. Rajagopal and U. A. Wiedemann, "Wilson loops in heavy 
ion collisions and thei r calculation in AdS/CFT," JHEP 0703 (2007) 066 
|arXiv:hep-ph/0612168] . 

[28] M. Chernicoff, J. A. Garcia and A. Giiijosa, "The energy of a moving 
quark-antiquark pair in an M = 4 SYM plasma," JHEP 0609 (2006) 068 
|arXiv:hep-th/0607089] . 

[29] P. C. Argyres, M. Edalati and J. F. Vazquez-Poritz, "Spacelike strings and jet 
quenching from a Wilson loop," JHEP 0704 (2007) 049 [arXiv:hep-th/0612157| . 

[30] P. C. Argyres, M. Edalati and J. F. Vazquez-Poritz, "Lightlike Wilson loops from 
AdS /CFT." rarXiv:0801.4594l [hep-th], 

[31] K. Peeters, J. Sonnenschein and M. Zamaklar, "Holographic melting and related 
properties of mesons in a quark gluon plasma," Phys. Rev. D 74 (2006) 106008 
|arXiv:hep-th/060619"5] . 

[32] H. Liu, K. Rajagopal and U. A. Wiedemann, "An AdS/CFT calcula- 
tion of screening in a hot wind," Phys. Rev. Lett. 98 (2007) 182301 
|arXiv:hep-ph/06070'62] . 

[33] T. Sakai and S. Sugimoto, "Low energy hadron physics in holographic QCD," 
Prog. Theor. Phys. 113 (2005) 843 |arXiv:hep-th/0412141| ; 
"More on a holograph ic dual of QCD," Prog. Theor. Phys. 114 (2006) 1083 
|arXiv:hep-th/050707"3] . 



59 



[34] E. Caceres, M. Natsuume and T. Okamura, "Screening length in plasma winds," 



JHEP 0610 (2006) Oil arXiv:hep-th/0607233|; 



S. D. Avramis, K. Sfetsos and K. Siampos, "Stability of strings dual to flux 
tubes between static quarks in N=4 SYM," Nucl. Phys. B 769 (2007) 44 
|arXiv:hep-th/0612139| ; 

C. Hoyos-Badajoz, K. Landsteiner and S. Montero, "Holographic Meson Melt- 



ing," JHEP 0704 (2007) 031 arXiv:hep-th/0612169| 



0. Antipin, P. Burikham and J. Li, "Effective Quark Antiquark Potential in 
the Quark Gluon Plasm a from Gravity Dual Models," JHEP 0706 (2007) 046 
|arXiv:hep-ph/0703105] ; 

M. Natsuume and T. Okamura, "Screening length and the direction of plasma 

winds," JHEP 0709 (2007) 039 |arXiv:0706.0086l [hep-th]]: 

S. D. Avramis, K. Sfetsos and K. Siampos, "Stability of string configurations dual 

to quarkonium states in AdS/CFT," Nucl. Phys. B 793 (2008) 1 [ arXiv:0706.2655l 

[hep-th]]; 

1. Amado, C. Hoyos-Badajoz, K. Landsteiner and S. Montero, "Absorption 
Lengths in the Holographic Plasma," JHEP 0709 (2007) 057 [arXiv: 0706 .27501 
[hep-th]]; 

H. Dorn and T. H. Ngo, "On the internal space dependence of the static quark- 
antiquark potential in JV = 4 SYM plasma wind," Phys. Lett. B 654 (2007) 41 
|arXiv:0707.2754l [hep-th]]; 

T. Song, Y. Park, S. H. Lee and C. Y. Wong, "The thermal width of heavy 
quarkonia moving in quark gluon plasma," Phys. Lett. B 659 (2008) 621 
|arXiv:0709.0794l [hep-ph]]; 

D. Hou and H. C. Ren, "Heavy Quarkonium States with the Holographic Poten- 
tial," JHEP 0801 (2008) 029 |arXiv:0710.2639l [hep-ph]]: 

K. Peeters and M. Zamaklar, "Dissociation by acceleration," JHEP 0801 (2008) 
038 |arXiv:0711.3446l [hep-th]]; 

E. Caceres, R. Flauger, M. Ihl and T. Wrase, "New Supergravity Backgrounds 
Dual to N=l SQCD-like Theories with N f = 2iV c ," larXiv: 071 1.48781 [hep-th]; 
K. D. Jensen, A. Karch and J. Price, "Strongly bound mesons at finite temper- 
ature and in magnetic fields from AdS/CFT," larXiv: 0801 .24011 [hep-th]. 

[35] P. C. Argyres, M. Edalati and J. F. Vazquez-Poritz, "No-drag string configura- 
tions for steadily moving quark-antiquark pairs in a thermal bath," JHEP 0701 



(2007) 105 |arXiv:hep-th/0608118|. 



[36] S. D. Avramis, K. Sfetsos and D. Zoakos, "On the velocity and chemical-potential 
dependence of the heavy-quark interaction in N = 4 SYM plasmas," Phys. Rev. 



D 75 (2007) 025009 arXiv:hep-th/0609079|. 



[37] J. J. Friess, S. S. Gubser, G. Michalogiorgakis and S. S. Pufu, "Stability of strings 



binding heavy-quark mesons," JHEP 0704 (2007) 079 arXiv:hep-th/0609137|. 



60 



[38] D. Bak, A. Karch and L. G. Yaffe, "Debye screening in strongly coupled N=4 
supersymmetric Yang-Mills plasma," JHEP 0708 (2007) 049 [ arXiv: 0705 .09941 
[hep-th]]. 

[39] Q. J. Ejaz, T. Faulkner, H. Liu, K. Rajagopal and U. A. Wiedemann, "A 
limiting velocity for quarkonium propagation in a strongly coupled plasma via 
AdS/CFT," larXiv:0712.0590l [hep-th]. 

[40] Y. Hatta, E. Iancu and A. H. Mueller, "Jet evolution in the N=4 SYM plasma 
at strong coupling," larXiv:0803.2481l [hep-th]. 

[41] C. G. Callan and A. Giiijosa, "Undulating strings and gauge theory waves," 
Nucl. Phys. B 565, 157 (2000) |arXiv:hep-th/9906153| . 

[42] S. Lin and E. Shuryak, "Stress Tensor of Static Dipoles in strongly coupled H—A 
Gauge Theory," Phys. Rev. D 76 (2007) 085014 jarXiv:0707.3T35l [hep-th]]. 

[43] S. Caron-Huot, P. Kovtun, G. D. Moore, A. Starinets and L. G. Yaffe, "Photon 
and dilepton production in supersymmetric Yang-Mills plasma," JHEP 0612 
(2006) 015 |arXiv:hep-th/0607237l ; 

A. Parnachev and D. A. Sahakyan, "Photoemission with chemical potential from 
QCD gravity dual," Nucl. Phys. B 768 (2007) 177 [arXiv:hep-th/0610247| ; 
D. Mateos and L. Patino, "Bright branes for strongly coupled plasmas," JHEP 
0711 (2007) 025 [arXiv: 0709 .21681 [hep-th]]; 

A. Nata Atmaja and K. Schalm, "Photon and Dilepton Production in Soft Wall 
AdS/QCD," larXiv:0802. 14601 [hep-th]. 

[44] Y. Hatta, E. Iancu and A. H. Mueller, "Deep inelastic scattering off a M = 4 
SYM plasma at strong coupling," JHEP 0801 (2008) 063 |arXiv:0710.5297l [hep- 
th]]. 

[45] S. S. Gubser, "Relativistic heavy ion collisions and string theory," lectures at 
Prospects in Theoretical Physics, 2006, available at 



http : //www . admin . ias . edu/pitp/2006f iles/Lecture/ 20notes/ 

Gubser°/ 20Lecture7o20notes7o20REVISED . pdf 

[46] S. S. Gubser, "Heavy ion collisions and black hole dynamics," Gen. Rel. Grav. 
39 (2007) 1533; 

M. Natsuume, "String theory and quark-gluon plasma," arXi v: hep-ph / 70 1 2 1 [ 
H. Liu, "Heavy ion co llisions and AdS/CFT," J. Phys. G 34 (2007) S361 
|arXiv:hep-ph/07022i0| ; 

K. Peeters and M. Zamaklar, "The string/gauge theory correspondence in QCD," 
larXiv: 0708 .15021 [hep-ph]; 

D. Mateos, "String Theory and Quantum Chromodynamics," Class. Quant. 
Grav. 24 (2007) S713 |arXiv:0709. 15231 [hep-th]]: 

J. Erdmenger, N. Evans, I. Kirsch and E. Threlfall, "Mesons in Gauge/Gravity 



61 



Duals - A Review," larXiv: 071 1.44671 [hep-th]; 

M. Natsuume, "String theory for heavy-ion physicists," Prog. Theor. Phys. 
Suppl. 168 (2007) 372. 

[47] D. Mateos, R. C. Myers and R. M. Thomson, "Thermodynamics of the brane," 
JHEP 0705 (2007) 067 [arXiv:hep-th/0701132] ; 

A. Buchel, S. Deakin, P. Kerner and J. T. Liu, "Thermodynamics of the N = 2* 
strongly coupled plasma," Nucl. Phys. B 784 (2007) 72 |arXiv:hep-th/0701142| ; 
O. Aharony, A. Buchel and P. Kerner, "The black hole in the throat - thermo- 
dynamics of strongly coupled cascading gauge theories," Phys. Rev. D 76 (2007) 
086005 [arXiv:0706.1768l [hep-th]]; 

M. Mahato, L. A. P. Zayas and C. A. Terrero-Escalante, "Black Holes in Cascad- 
ing Theories: Confinement /Deconfmement Transition and other Thermal Prop- 
erties," JHEP 0709 (2007) 083 |arXiv:0707.2737l [hep-th]]; 
J. Schmude, "The quark-gluon plasma and D6-branes on the conifold," 
larXiv:0711.3763l [hep-th]. 

[48] S. S. Gubser, "Comparing the drag force on heavy quarks in M = 4 super- Yang- 
Mills theory and QCD," Phys. Rev. D 76 (2007) 126003 |arXiv:hep-th/0611272|. 



[49] H. Nastase, "The RHIC fireball as a dual black hole," |arXiv:hep-th/0501068 



E. Shuryak, S. J. Sin and I. Zahed, "A gravity dual of RHIC collisions," J. Ko- 
rean Phys. Soc. 50, 384 (2007) |arXiv:hep-th/0511199] ; 

R. A. Janik and R. B. Peschanski, "Asymptotic perfect fluid dynamics as a con- 
sequence of AdS/CFT," Phys. Rev. D 73, 045013 (2006) |arXiv:hep-th/0512162| ; 
H. Nastase, "More on the RHIC fireball and dual black holes," 



arXiv:hep-th/0603176 



R. A. Janik and R. B. Peschanski, "Gauge / gravity duality and thermal- 
ization of a boost-invariant perfect fluid," Phys. Rev. D 74, 046007 (2006) 



arXiv:hep-th/0606149|; 



S. Nakamura and S. J. Sin, "A holographic dual of hydrodynamics," JHEP 0609, 
020 (2006) |arXiv:hep-th/0607123| ; 

S. J. Sin, S. Nakamura and S. P. Kim, "Elliptic flow, Kasner universe and holo- 
graphic dual of RHIC fireball," JHEP 0612, 075 (2006) |arXiv:hep-th/0610113| ; 
R. A. Janik, "Viscous plasma evolution from gravity using AdS/CFT," Phys. 
Rev. Lett. 98, 022302 (2007) |arXiv:hep-th/0610144| ; 

S. Lin and E. Shuryak, "Toward the AdS/CFT gravity dual for High Energy 



Collisions: I.Falling into the AdS," |arXiv:hep-ph/0610168t 



J. J. Friess, S. S. Gubser, G. Michalogiorgakis and S. S. Pufu, "Expanding plas- 
mas and quasinormal modes of anti-de Sitter black holes," JHEP 0704 (2007) 
080 |arXiv:hep-th/0611005| ; 

D. Bak and R. A. Janik, "From static to evolving geometries: R- 
charged hydrodynamics from supergravity," Phys. Lett. B 645 (2007) 303 
|arXiv:hep-th/0611304| ; 

K. Kajantie and T. Tahkokallio, "Spherically expanding matter in AdS/CFT," 



62 



Phys. Rev. D 75 (2007) 066003 [arXiv:hep-th/0612226| ; 



M. P. Heller and R. A. Janik, "Viscous hydrodynamics relaxation time from 
AdS/CFT," Phys. Rev. D 76 (2007) 025027 |arXiv:hep-th/0703243| ; 
K. Y. Kim, S. J. Sin and I. Zahed, "Diffusion in an Expanding Plasma using 
AdS/CFT," llarXiv:0707.060Tl [hep-th]; 

D. Bak, M. Gutperle and A. Karch, "Time dependent black holes and thermal 
equilibration," JHEP 0712 (2007) 034 |arXiv:0708.3691l [hep-th]]: 
J. Grosse, R. A. Janik and P. Surowka, "Flavors in an expanding plasma," 
larXiv:0709.39T0l [hep-th]; 

S. Lin and E. Shuryak, "Toward the AdS/CFT Gravity Dual for High Energy 
Collisions: II. The Stress Tensor on the Boundary," larXiv:0711.0 736 [hep-th]; 
P. Benincasa, A. Buchel, M. P. Heller and R. A. Janik, "On the supergravity 
description of boost invariant conformal plasma at strong coupling," Phys. Rev. 
D 77 (2008) 046006 |arXiv:0712.2"025l [hep-th]]: 

J. Alsup and G. Siopsis, "Bjorken flow from an AdS Schwarzschild black hole," 
larXiv:0712.2To4l [hep-th]; 

A. J. Amsel, D. Marolf and A. Virmani, "Collisions with Black Holes and De- 
confined Plasmas," larXiv:0712.222Tl [hep-th]; 

S. Bhattacharyya, V. E. Hubeny, S. Minwalla and M. Rangamani, "Nonlinear 
Fluid Dynamics from Gravity," JHEP 0802 (2008) 045 |arXiv:0712.2"456l [hep- 
th]]; 

M. Natsuume and T. Okamura, "Comment on "Viscous hydrodynamics relax- 
ation time from AdS/CFT correspondence'," larXiv:0712.2 917 [hep-th]. 

[50] S. Peigne, P. B. Gossiaux and T. Gousset, "Retardation effect for collisional 
energy loss of hard partons produced in a QGP," JHEP 0604 (2006) 011 
|arXiv:hep-ph/0509185] . 

[51] M. Djordjevic, "Collisional energy loss in a finite size QCD matter," Phys. Rev. 
C 74 (2006) 064907 |arXiv:nucl-th/0603066| ; 

P. B. Gossiaux, S. Peigne, C. Brandt and J. Aichelin, "Energy Loss of a 
Heavy Quark Produce d in a Finite Size Medium," JHEP 0704 (2007) 012 
|arXiv:hep-ph/06080'6l] ; 

P. B. Gossiaux, J. Aichelin, C. Brandt, T. Gousset and S. Peigne, "Energy loss 
of a heavy quark produced in a finite-size quark-gluon plasma," J. Phys. G 34 
(2007) S817 |arXiv:hep-ph/0703095| . 

[52] A. Mikhailov, "Nonlinear waves in AdS/CFT correspondence," 



|arXiv:hep-th/0305196 



[53] A. Karch and E. Katz , "Adding flavor to AdS/CFT," JHEP 0206 (2002) 043 
|arXiv:hep-th/020523"6] . 

[54] H. Arodz, "Solutions Of Yang-Mills Equations With External Sources," Phys. 
Lett. B 78 (1978) 129; 



63 



A. Trautman, "Radiation Of Energy And Change In Color Of A Point Source 
Of The Yang-Mills Field," Phys. Rev. Lett. 46 (1981) 875; 
W. Drechsler, "Does A Moving Nonabelian Charge Radiate?," Phys. Lett. B 90, 
258 (1980); 

W. Drechsler and A. Rosenblum, "Equations Of Motion And Iteration Of 
Lienard-Wiechert Type Solutions In Classical Yang-Mills Theory," Phys. Lett. 
B 106 (1981) 81; 

S. Schlieder, "Some Remarks On Charges And Their Conservation In A Classical 
SU(2) Yang-Mills Theory," Nuovo Cim. A 63, 137 (1981); 

J. Tafel and A. Trautman, "Can Poles Change Color?," J. Math. Phys. 24, 1087 
(1983); 

C. H. Lai and C. H. Oh, "Color Screening In Classical Yang-Mills Theories With 
Sources," Phys. Rev. D 29, 1805 (1984); 

C. H. Oh, C. H. Lai and R. Teh, "Color Radiation In The Classical Yang-Mills 
Theory," Phys. Rev. D 33 (1986) 1133; 

O. Sarioglu, "Lienard-Wiechert potentials of a non-Abelian Yang-Mills charge," 
Phys. Rev. D 66 (2002) 085005 |arXiv:hep-th/0207227| ; 

L. D. McLerran and R. Venugopalan, "Computing quark and gluon dis- 
tribution functions for very large nuclei," Phys. Rev. D 49 (1994) 2233 
|arXiv:hep-ph/9309289| . 

[55] J. F. Gunion and G. Bertsch, "Hadronization By Color Bremsstrahlung," Phys. 
Rev. D 25 (1982) 746; 

R. D. Field, G. C. Fox and R. L. Kelly, "Gluon Bremsstrahlung Effects In Hadron 
- Hadron Collisions," Phys. Lett. B 119 (1982) 439; 

M. Gyulassy and L. D. McLerran, "Yang-Mills radiation in ultrarelativistic nu- 
clear collisions," Phys. Rev. C 56 (1997) 2219 |arXiv:nucl-th/9704034| . 

[56] U. H. Danielsson, E. Keski-Vakkuri and M. Kruczenski, "Vacua, propa- 
gators, and holograp hic probes in AdS/CFT," JHEP 9901, 002 (1999) 
|arXiv:hep-th/981200"7| . 

[57] K. Maeda and T. Okamura, "Radiation from an accelerated quark in AdS/CFT," 
larXiv: 0712 .41201 [hep-th]. 

[58] L. Susskind and E. Witten, "The Holographic Bound In Anti-De Sitter Space," 



arXiv:hep-th/9805114 



A. W. Peet and J. Polchinski, "UV/IR relations in AdS dynamics," Phys. Rev. 
D 59 (1999) 065011 |arXiv:hep-th/9809022| . 

[59] J. L. Hovdebo, M. Kruczenski, D. Mateos, R. C. Myers and D. J. Winters, 
"Holographic mesons: Adding flavor to the AdS/CFT duality," Int. J. Mod. 
Phys. A 20 (2005) 3428. 



64 



[60] S. Kobayashi, D. Mateos, S. Matsuura, R. C. Myers and R. M. Thomson, 
"Holographic phase transitions at finite baryon density," JHEP 0702 (2007) 
016 |arXiv:hep-th/0611099| . 

[61] L. Brink and M. Henneaux, "Principles of String Theory, 11 New York, USA: 
Plenum (1988) 297p 

[62] J. Gomis and F. Passer ini, "Holographic Wilson loops," JHEP 0608, 074 (2006) 
|arXiv:hep-th/0604007| . 

[63] S. J. Sin and I. Zahed, "Holography of radiation and jet quenching," Phys. Lett. 
B 608 (2005) 265 |arXiv:hep-th/0407215| . 

[64] S. J. Rey, S. Theisen and J. T. Yee, "Wilson-Polyakov loop at finite temperature 
in large N gauge theory and anti-de Sitter super gravity," Nucl. Phys. B 527 
(1998) 171 |arXiv:hep-th/9803135l ; 

A. Brandhuber, N. Itzhaki, J. Sonnenschein and S. Yankielowicz, "Wilson loops 
in the large N limit at finite temperature," Phys. Lett. B 434, 36 (1998) 
|arXiv:hep-th/980313"7| . 

[65] J. M. Maldacena, "Wilson loops in large N field theories," Phys. Rev. Lett. 80 
(1998) 4859 [arXiv:hep-th/9803002] ; 

S. J. Rey and J. T. Yee, "Macroscopic strings as heavy quarks in large N 
gauge theory and anti-de Sitter super gravity," Eur. Phys. J. C 22 (2001) 379 
|arXiv:hep-th/980300l] . 

[66] M. Laine, O. Philipsen, P. Romatschke and M. Tassler, "Real-time static poten- 



tial in hot QCD," JHEP 0703 (2007) 054 arXiv:hep-ph/0611300 



[67] G. D. Moore and D. Teaney, "How much do heavy quarks thermalize in a heavy 
ion collision?," Phys. Rev. C 71 (2005) 064904 |arXiv:hep-ph/0412346] . 

[68] I. R. Klebanov, W. Taylor and M. Van Raamsdonk, "Absorption of dilaton par- 
tial waves by D3-branes," Nucl. Phys. B 560 (1999) 207 [arXiv:hep-th/9905174| . 

[69] M. Kruczenski, D. Mateos, R. C. Myers and D. J. Winters, "Meson spectroscopy 
in AdS/CFT with flavour," JHEP 0307 (2003) 049 [arXiv:hep-th/0304032] . 

[70] S. Hong, S. Yoon and M. J. Strassler, "Quarkonium from the fifth dimension," 
JHEP 0404 (2004) 046 [arXiv:hep-th/0312071] . 

[71] H. Satz, "Quarkonium binding and dissociation: The spectral analysis of the 
QGP," Nucl. Phys. A 783 (2007) 249 |arXiv:hep-ph/0609197| . 

[72] O. Philipsen, "Non-perturbative formulation of the static color octet potential," 
Phys. Lett. B 535, 138 (2002) |arXiv:hep-lat/0203018| ; 

O. Kaczmarek, F. Karsch, P. Petreczky and F. Zantow, "Heavy quark anti-quark 



65 



free energy and the renormalized Polyakov loop," Phys. Lett. B 543, 41 (2002) 
|arXiv:hep-lat /0207002| . 



6(3 



